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PART  I     . 

.-.LYSIS  AND   OPTIMIZATION   OF  STEP  AERATION 
ACTIVATED   SLUDGE   PROCESS 


CHAPTER  I 
INTRODUCTION 
I.    BIOCHEMICAL  ENGINEERING  AND  CHEMICAL  ENGINEL.  .. 

Various  definitions  of  biochemical  engineering  have  been 
given.   H.  Hartley  (1)*  has  said  that  "Biochemical  engineer- 
ing is  the  practical  application  of  our  knowledge  of  micro- 
nanism.   The  industrial  application  of  a  biochemical 
orocess  is  a  combined  operation  in  which  the  biochemist, 
.  b  microbiologist,  the  geneticist,  and  the  chemical  en- 
gineer are  all  intimately  concerned."   In  lectures  on  bio- 
chemical engineering  at  the  University  of  Tokyo  in  1963, 
the  following  definition  has  been  offered.   "Biochemical 
engineering  is  that  activity  concerned  with  economic  pro- 
cessing of  materials  of  biological  character  or  origin  to 
serve  useful  purposes.   The  function  of  the  biochemical 
engineer  is  that  of  translating  the  knowledge  of  the  mi- 
crobiologist and  Lhe  biochemist  into  a  practical  operation." 

From  the  foregoing  definition,  it  is  obvious  that  bio- 
chemical engineering  is  intimately  related  to  chemical 

rir.eering.   A  biochemical  engineer  must  not  only  have  an 
appreciation  of  the  biological  science,  but  he  also  must 
be  well  grounded  in  chemical  engineering  principles. 


—T.'cers  in  parentheses  refer  to  references  given  on 
page  143. 


ATER  T. ENT  AND  ITS  PROCESSES 

tnicipal  and  industrial  wastes  are  treated  to  protect 
pu       ilth,  to  avoid  nuisance,  to  prevent  the  polluting 
of  na       iter  arc  of  industrial  water,  and  to  avoid  dam- 
5  suits.   The  pollutional  characteristics  of  waste  waters 
e  classified  according  to  their  state  (suspended,  ccl- 
L,  and  dissolved)  and  their  nature  (inorganic,  organic, 
gases,  and  living  organisms).   However  the  major  components 
of  organic  wastes  are  usually  the  suspended  solids  and  the 
organic  content.   The  treatment  of  this  waste  is  usually 
carried  out  in  four  steps,  viz.  pretreatment ,  biological 
oxioa'icr. ,  sludge  treatment,  and  disposal. 

Pretreatment  includes  screening,  grit  removal,  arc 
sedimentation  or  flotation.   Biological  oxidation  is  usually 
accomplished  in  either  fixed-bed  units  (trickling  filters) 
or  in  fluid-bed  systems  (activated  sludge).   It  is  employed 
in  removing  the  colloidal  and  dissolved  organic  matter  and 
it  also  olays  the  most  important  role  in  waste  treatment. 

j  from  the  foregoing  units  requires  further  treatment, 
on  and  dewatering  before  disposal.   Let  us  consider 
»te  treatment  plant  as  shown  in  Figure  1.   Generally, 
:ed  sludge  process  is  a  part  of  such  a  plant.   The 
Ling  tank  removes  settleable  solids  from  the  raw 
stre^  .   The  settled  solids  are  diverted  for  stabili- 
se anaerobic  digestion.   The  waste  stream  from 
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settlii    .  :ers  the  aeration  phase  where  the  or- 
.  .  food  source  for  aerobic  micro- 
anisms  convert  a  part  02"  the  total  food 
:c  biological  cell  mass  (this  is  the  synthesis  step) 
Lch  remains  in  the  liquid  stream  leaving  the  aeration  tank 
fraction  of  the  total  food  is  utilized  by  the 
organisms  as  a  source  of  energy  for  cell  synthesis  and 

:e.   The  latter  fraction  of  the  food  is  converted 
to  carbon  dioxide  and  water.   The  settleable  biological 
sclids     Lnerts  are  removed  from  the  waste  stream  in  the 
final  settling  tank.   A  part  of  this  activated  sludge  is 
returned  to  the  aeration  tank  to  maintain  a  favorable  bio- 
_   ical  solids  concentration  in  it.   The  excess  solids 
st  be  stabilized  further,  often  by  anaerobic  digestion. 

3.    ACTIVATED  SLUDGE  PROCESS 

The  activated  sludge  process  may  be  defined  as  a  system 
in  which  flocculated  biological  growths  are  continuously 
circulated  and  contracted  with  organic  waste  in  the  presence 
of  oxygen.   The  oxygen  is  usually  supplied  from  air  bubbles 
ected  into  the  sludge-liquid  mass.   The  process  involves 
aeration  step  followed  by  a  solid-liquid  separation  step 
h  the  separated  sludge  is  recycled  back  and  mixed 
waste.   A  portion  of  this  sludge  is  re- 
.  .  ....  ..-.  1  fcreatm  nt  and  disposal.   The  aeration  step 

•  be  considered  in  three  functional  phases  as  follows  (2): 


(a)  A  rapid  adsorption  of  waste  substrate  by  the 
activated  sludge. 

(b)  Progressive  oxidation  and  synthesis  of  the  ab- 
sorbed organics. 

(c)  Further  aeration  resulting  in  oxidation  and  dis- 
persion of  the  sludge  particles. 

Various  modifications  of  the  activated  sludge  process 
have  been  developed  to  achieve  economic  advantage  in  con- 
struction and  operation  (3).   The  conventional  activated 
sludge  process  as  shown  in  Figure  2  is  one  of  the  basic 
activated  sludge  processes.   The  present  knowledge  con- 
cerning the  treatment  has  been  summarized  very  effectively 
by  Kaseltine  (4).   According  to  his  report,  there  are  some 
limitations  to  this  process.   In  order  to  overcome  the  dis- 
advantages of  the  process,  some  modifications  of  the  conven- 
tional process  have  been  proposed.   The  step  aeration  system 
as  shown  in  Figure  3,  developed  by  Gould  (5),  is  one  of  the 
modifications.   In  essence,  the  modified  process  involves 
the  feeding  of  the  sewage  flow  at  different  points  along 
the  aeration  basin.   The  organic  load  is  thus  distributed 
over  the  length  of  the  basin.  As  a  result,  accelerated 
growth  and  oxidation  are  not  confined  to  one  end  of  the  basin 
as  in  the  conventional  process  but,  instead,  take  place  over 
most  of  the  basin.   The  difference  between  the  oxygen  demand 
patterns  of  the  conventional  aeration  and  step  aeration  is 
illustrated  in  Figures  3  and  l+. 
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A.        IGN  AND  .       JS 

Economy  has  always  been  an  important  engineering  con- 
sideration in  any  process  industry.   A  successful  process 
design  will  permit  greater  economy  in  yielding  products. 
Optimization  in  turn  plays  a  very  important  role  in  process 
desir-r. .   In  the  study  of  the  step  aeration  process,  the 
distribution  of  the  organic  load  and  the  flow  behavior  are 
found  to  be  important  variables  which  affect  the  optimal 
design  of  the  process.   In  this  study,  four  types  of  systems 
are  investigated  for  the  overall  pattern  of  flow,  mixing, 
and  distribution  of  feed  for  the  step  aeration  process. 
The  combination  of  the  discrete  version  of  the  maximum 
principle  and  the  steepest  descent  search  technique  is  used 
to  carry  out  the  optimization  work. 


lo 


chapter  ii 
lTICAL  analysis  of  the  process 

AT]  D  KINETIC  MODEL 

(a)    rinciple  of  Biological  Oxidation  (6) 

oiogical  oxidation  is  simply  a  conversion  process 
therein  dissolved  organic  compounds  are  converted  into  bac- 
terial cells,  which  can  then  be  removed  from  the  waste 

or.   The  generalized  reaction  for  the  removal  of  soluble 
organics  may  be  considered  as  follows: 


.  -ar.ics  +  02  +  Nutrients  Organisms  >  CG^   +  H  0  +  Organisms 


(b)   Growth  Pattern  (2) 

3  curve  in  Figure  5  illustrates  the  classic  growth 
item  exhibited  by  microorganisms  in  a  batch  culture.   Ex- 
tion  of   the  curve  reveals  that  growth  passes  through  ■ 
»rent  phases.   Initially,  all  nutrients  are  pre- 
sent in  excess  of  the  requirements  of  the  microorganism, 
.  growth  is  unrestricted.   During  this  period,  called  the 
,ant  growth  phase,  the  concentration  of  microorganisms 
incj  -  .s  at  an  exponential  rate.   At  some  concentration, 

nutrients  becomes  growth  limiting  and  the  culture 
proceeds  into  the  declining  growth  phase.   In  response  to 

competition  of  the  microorganisms  for  the 
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Fig.  5.    Classic    growth   pattern 
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nutrient,  the  rate  of  growth  decreases 
Lily  halts.   The  remaining  portion  of  the 
the  decrease  of  the  microorganism  resulting 
autooxidation  which  occurs  after  the  depletion  of  the 

nics.   This  is  often  called  the  endogeneous 

res-      ..  phase  of  activated  sludge. 

(c)   Kinetic  Model 

As  stated  previously,  in  an  aerobic  environment  net 
limited  by  mixing  or  nutrients,  the  specific  growth  rate 
is  proportional  to  the  concentration  of  organisms.   The 
mat]    ;ical  model,  therefore,  can  be  represented  by 


QX/, 


KXj 


dt      ~%  (1) 

,re 

r  =  growth  rate,  mg/liter,  hr., 
Xi  =  concentration  of  organisms,  mg/liter, 
k  =  growth  rate  constant,  hr"  . 

nod  (7),  in  his  study  of  bacterial  growth  in  a  chemostat 
(cr     iouj  culture),  found  that  his  data  behaved  in  a  man- 
Led  by  the  Michaelis-Menten  equation  which  is  one 
of  t;     st  widely  accepted  models  of  enzyme  kinetics.   He 
re  su     :ed  the  following  form  for  the  growth  rate 
c  o  i    n  t . 
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k  =  k  s   (   x2   }  (2) 

...ax  K  -1-  ^2 


3re 

k    =  maximum  growth  rate  when  the  organic  concen- 

_~x  " 

t rat ion  is  not  limiting  the  rate  of  growth,  hr~  , 
K  =  the  concentration  of  crganics  at  which  the 

specific  growth  rate  observed  is  one  half  the 
maximum  value;  saturation  constant,  mg/liter, 
Xp  =•  concentration  of  organics,  mg/liter. 

Then,  if  one  milligram  of  organics  produces  Y  milligrams  of 
organisms,  the  rate  of  consumption  of  organics  will  be 

s   dt      Y  dt 

From  an  engineering  viewpoint,  we  may  consider  the 
various  phases  of  sludge  growth  and  BOD  removal  to  consist 
of  a  dynamic  relationship  between  the  mass  transfer  of  es- 
sential foods  into  the  cell  structure,  and  the  assimilation 
and  utilization  of  these  foods  for  energy  and  growth.   At 

jh  concentrations  of  organic  matter,  the  rate  of  assimi- 
lation and  the  growth  rate  are  independent  of  the  concen- 
tration of  organic  matter.   At  low  organic  levels,  the  rate 
of  growth  and  hence  the  rate  of  BCD  removal  are  frequently 
concentration  dependent.   The  kinetic  form  of  Equation  (2) 
describes  the  situation  except  during  the  endogenous  phase. 
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r  the     ~  phase  is  independent  of  the 
ion  of  organics,  X£,  since  the  organics  are  in 
excess  i.e.  X£>>  K,  and  thus  Equation  (2)  becomes 

.  .^x 

During  the  initial  part  of  the  second  phase  of  growth,  k 
.  on  X2  and  the  original  form  of  Equation  (2) 
plies,  but  as  the  organics  are  consumed,  x^  becomes  very 
ill,  i.e.  X£  <<  K  and  Equation  (2)  simplifies  to 

=  k;rnx  *2  .  U) 

.. 

During  the  endogeneous  phase,  the  effect  of  endogenecus 
respiration  becomes  important.   In  the  expression  for  one 
phase  kinetics  proposed  by  Grieves  et  al.  (8),  the  effect 
endogeneous  respiration  is  considered.   The  expression 
is 


-  -'7'ax   *2  x,  -  kDx,  (5) 

x2  ^ 


re 

e  specific  enaogeneous  microbial  attrition  rate, 
-] 

,  equation  (5)  becomes 
■    .  ...  K  -k°x4-  (lc*-  kD)x^  ,  (6) 
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.  when  x„  «    K,  Equation  (5)  becomes 


— 
K  ^ 


k~  QV  XqX  i      n  . 

r  =  • ,\  -  -  -  kux,  .  (7) 


Theoretically,  when 


Xo=A 


max 


we  have 

r  =  0. 

Numerically  yp   is  usually  two  orders  of  magnitude  less  than 
^rrf  .   Hence  the  effect  of  endogeneous  respiration  is  neg- 

I* id  .A. 

legible  except  when  x?  is  very  small.   Therefore,  the  model 
represented  by  Equation  (2)  is  used  in  the  present  work. 

2.    THE  MATHEMATICAL  REPRESENTATION  OF  THE  PROCESS 
(a )   Flow  Models 

A  plug  flew  model  is  often  used  to  represent  the  flow 
behavior  of  the  conventional  activated  sludge  process.   This 
model  assumes  that  there  is  no  longitudinal  mixing  in  the 
aeration  tank.   However,  the  optimum  degree  of  longitudinal 
mixing  in  the  aeration  tank  is  still  in  question.   Some  in- 
vestigator (9)  has  advocated  the  use  of  aeration  tanks  de- 
signed to  provide  a  uniform  composition  within  the  entire 
tar..  .  Tr.e   latter  system  is  often  referred  to  as  the  complete 
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activated  sludge  process  (9)  and  the  completely 

.  .  stirred  tank  model  is  usually  used  to  represent  the 
f]       rlor   in  this  system. 

In  step  aeration,  the  influent  is  fed  to  the  reactor 
:  .     several  different  points  in  order  to  avoid  the 
confinement  of  growth  and  oxidation  to  one  end  of  the  basin 
he  conventional  process.   Furthermore,  from  the 
tiic  model  we  can  see  that  feeding  a  larger  amount  of 
feed   to  the  front  end  will  increase  the  organic  concen- 
tration but  dilute  the  organism  concentration.   Therefore, 
a  step  aeration  process  with  appropriate  distribution  of 
organic  load  may  perform  better  than  the  conventional  pro- 
cess.  A  series  of  plug  flow  units  in  which  feed  enters  at 
the  inlet  of  each  plug  flow  unit  can  be  used  to  represent 
■  flow  behavior  of  the  process.   Similarly,  the  step 
ration  process  can  also  be  represented  by  a  sequence  of 
continuously  stirred  tanks  in  which  feed  enters  the  system 
at  several  points. 

tplifying  Assumptions 

.'ollowing  assumptions  and  simplifications  are  made 

.  :  the  process  and  developing  the  mathematical 
ttion  for  the  process. 
(1)    te  system  is  isothermal  and  is  under  the  steady 
state  condition. 
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(2)  Physical  properties  such  as  density,  diffusivity, 
ar.ci  viscosity  are  constant. 

(3)  Y  is  dependent  only  upon  the  property  of  the  wacte 
itself,  and  independent  of  the  age  of  the  organ- 
isms and  the  effect  of  other  physical  conditions 
such  as  the  concentrations  of  organics  and  or- 
ganisms. 

(l+)      Organics  and  organisms  are  distinctly  separate 
entities  in  solution. 

(5)  Endogeneous  respiration  does  not  influence  the 
systerr.  performance. 

(6)  The  sludge  and  waste  streams  are  completely  mixed 
at  each  point  where  the  waste  is  introduced. 

(7)  Sufficient  oxygen  is  supplied  for  the  oxidation. 
(3)   The  fluid  is  a  continuum  and  there  is  no  segre- 
gation. 

Some  of  these  assumptions  depart  from  reality.   They 
are  justified  on  the  grounds  that  they  simplify  the  rela- 
tionship of  the  process  without  appreciably  changing  the 
"casic  characteristics  of  itself. 

(c)  r^'r.o   y.athematical  Representation 

sed  on  the  preceding  assumptions  and  the  flow  models 
described  in  the  previous  section  we  can  visualize  the 
system  as  consisting  of  mixing  stages  and  reaction  stages 
as  showr.  in  figure  6.   For  convenience,  we  consider  these 
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orer.t  parts  separately  ar.j  number  each  stage  u~ 
.gure  7-     also  define  the  following  notations 
for  the  subsequent  analysis. 

■  volumetric  rate  of  the  flow  discharged  fro:.,  the 
2n  th  stage,  liter/hr, 

x    =  concentration  of  organics  in  the  cutlet  stream 
2 

of  the  2n  th  stage,  mg/liter, 

X2   =  concentration  of  organics  in  feed,  mg/liter, 

2r 
Xo  *  ■»  accumulated  volume  up  to  and  including  the  2n  th 

stage,  liters, 

2r 
x,    =  concentration  of  organisms  in  the  outlet  stream 
4 


of  the  2n  th  stage,  mg/liter 


j 


x,        -   concentration  of  organisms  in  feed,  mg/liter, 

6-    =  volumetric  rate  of  waste  introduced  at  the 
(2n-l)th  stage,  liter/hr, 

- 

Q 2^     ~   volume  of  the  2n  th  stage,  liters. 

(1)   Analysis  of  a  mixing  stage 

In  order  to  establish  a  mathematical  model  of  the  mixing 
stage,  we  visualize  the  (2n-l)th  stage  as  shown  in  Figure 
6(a),  and  the  essential  governing  equations  are  written 
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Fig. 8     A    schematic    representation  of 
a  mixing  stage  and  a  reaction  stage, 
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he  assum    Lnetic  ar.c  flow  models  using  the  no- 
stablished  above, 
overall  material  baiar.c^. 


2n-2  ,  , 

x    =  y,    +  x,  (oj 

1      1 


._.   material  balance  of  the  oxidizable  orranics. 


2n-l  f    2n-2  2n-2    2n-l  2n-l 
01         x2  +  x2    x±         =x2    xx    .  (9) 

Equations  (6)  and  (9)  to  eliminate  xjr   and  solving 

for  the  exit  concentration  of  or^anics,  x£n~^,  we  obtain 

2n-l  f    2n-2  2n-2 

x2n-l  =  Gl         x2   *  x2         Xl  (10) 

2        g— --1  +  v2n-2 

6i      Xl 


.  organism  material  balance. 


..  2n-l       2n-2  2n-2    2n-l  2n-l  .  ■ 

o ,    x,  +  x,    x,    =  x,    x.,     .  (11) 

14    4     1       4     1 


is  assumed  to  be  equal  to  zero  and,  as  given  by 
Lon  (8) , 


2r.-.    2n-l    2n-2 
x     =  0     +  x    , 

1       1       1 


-1)  becoi 


.   .     2n-2  v2n-2 
•—  =  "I  (12) 

.2n-l    2n-2 


.v.  accumulated  volume 


x2n-l  =  x2n-2#  (13) 

3      3 


The  mixing  stage  is  assumed  to  mix  without  capacity. 

(2)  Analysis  of  a  reaction  stage 

A  reactor  stage  designated  as  the  2n  th  stage  is  shown 
in  Figure  8(b).   The  essential  governing  equations  are  de- 
rived below. 
i.   overall  material  balance. 


2n    2n-l  ,lM 

X]_  =  x-|_  (14  J 


ii.   organism  material  balance.   The  quantity  of  organisms 
produced  in  each  reactor  can  be  expressed  in  terms  of  the 
quantity  of  organics  consumed,  that  is, 

2n    2n-l  2n-l   ^  2n-l  2n-l    2n  2n 
Xl  XU     "  Xl    X4    "  Y(X1    X2    '  Xl  X2    }'   (15) 

ire  Y  denotes  the  yield  factor  of  the  biochemical  oxidation 
process.   Using  Equation  (14),  Equation  (15)  can  be  simplified 
to 
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-  :,'"*  ~  -  Y(x2  '    -  x2  )  .  (16) 


2n-l   ...    2n-l    2n 


.  >_ 


ecu.        volume. 


2n-l   Q2n  /.„, 

x~  =  x-    +  0O   .  (17) 


->     ^ 


2 


iv.  c.  anic  material  balance.  The  rate  of  growth  is  deter- 
.-.::.  the  concentration  of  organisms  and  organics  at  each 
point  within  the  reactor,  which  is  in  turn  influenced  by  the 

flow  model  used.   Therefore,  the  complete  mixing  model  and 
the  plug  flow  model  yield  two  different  material  balances 
for  the  oxidizable  organics. 

Since  the  exit  stream  from  the  continuous  flow  stirred 
tank  reaceor  has  the  same  concentration  of  organics  and  cr- 
nisms  as  the  fluid  within  the  reactor,  Equations  (1),  (2) 
and  (3)  car.  be  combined  to  give,  if  the  completely  stirred 
tank  system  is  considered, 

:2n    k  S      x|n      . 

-rs  -  -3 S^c  ( 2   j   2n  g  (18) 

Y   K  +  x2n   4 

ending  organic  material  balance  is 


-xf^-xf  xf  +  raef-0.  (19) 


iute  Equation  (io)  into  Equation  (19),  the  or- 


2$ 


,  s   2n  2n  2n 
-1  2n-l    2n   2n        x:-  x:  *  2  ,  nl 

Y(X  +  x2' 

For  this  system,  the  concentration  of  organisms,  x,  , 

is  dependent  upon  that  of  organics,  hence  xj"  can  be  expressed 

4 

as  function  of  x1.   By  applying  the  initial  stage  up  to  and 

including  the  ith  stage,  we  obtain 

Yx-j_  X2  -  Yx-j_  x^>  -  Yx£  x^ 


i          1 
x,    = 

4     xl 

CO        _  0 

0 
x2 

*1 

Yx^   -   Yx^ 
1X2         ^2    ' 

i  =  1,  2,  ...  2N,  (21) 


where 


A   =  x°  x°  +  Yx°  x°  -  Yx«  x^  '22) 

o     1  U  12      12 


Since  the  plug  flow  model  assumes  that  there  is  no 
mixing  of  fluid  longitudinally  along  the  flow  path,  the  con- 
centration of  organics  and  organisms  varies  from  position  to 
position  along  the  reactor.   Therefore,  the  material  balance 
for  organics  must  be  made  for  a  differential  element  of 
volume  d  . 

if erring  to  Figure  9,  the  material  balance  is 

(2n),,.   (2n)/0.     (2n)  ,  ,   ,  fl»    (2n)  ,  „   jnx     (2n)„  ,„ 
Xj_    \i)   x2    (4)  -  x1    (i+  dx)     x-?    (  -  +  o^ )  f  r<;   A  dJE 

-  0  (23) 
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where  x(2ri^(J2)  and  xj2n)(l  +  di)  are  the  volumetric  flow 

ra^es  aL  positions  j2  and  ^+  dJZ  respectively  of  the  2n  ... 

stage,  x    (£)  and  x^  n'(  /  +  d/)  are  the  concentration  of 

2 

.cs  at  positions  2.   and  j? -i-  di  respectively  of  the  2n  th 
s"cage,  and  r     is  the  reaction  rate  within  the  differential 
volume,  dV,  between  j?  and  Jl  +  dJL.      Since  density  is  constant, 
the  volumetric  flow  rate  is  unchanged  along  the  plug  flow 
reactor.   Therefore,  Equation  (23)  can  be  simplified  and 
beccrr.es 

-x(?n)(i)  dx<2rJ(£)  +  rl2n)Adi=  0  .  (24) 


Rearranging  this  equation  gives 


d42n^>   _         A  r(2n) 

xj2n,(i)    's 


(25) 


where 


-kx,(2n)(4)        -k  L  x^2n)(i)   x,(2n)(i) 


-Kr.ax  x2  ^   XL 

Y  (K  +  x^2nj(Jl)j 


r£    "     =  Y =  ~~"77T~ (2n)        ^ ■ 


f  2n ) 
re  x.v   '(J?)  is  the  concentration  of  organisms  at  the 

position  x  of  the  2n  th  stage,  which  is  dependent  upon 

x2~'         )•   Since  x|2n)  (jg )  is  constant  throughout  the  stage, 

Equation  (25)  can  be  integrated  if  x/2n^(i)  can  be 

4 

expressed  in  terms  of  x^2n'(i). 
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lire  10,  the  overall  orgar.is:.-.  material 
2€  - -'     te  1st  stage  up  zc   the  position  X  of  the  2r.  th 
be  use^  to  express  x*2n)(j2)  as 

"  (i)  =  a  -  Yxi2n)U),  (27) 


2n- 

\-  ~  -  • 


i     (xOxo  +  Yxoxo.yxoxf]  _ 


If  v;e  substitute  Equation  (27)  into  Equation  (25)  and  inte- 

-  applying  the  boundary  conditions, 


(2n)/rn    2n-i     .    (2n),T,    2n  ,„_. 

x-„    (vj)  =  x2     ana  x2    (iJ  =  x2  U9) 


Dbtian 


(a  -  Yx^"1)    v    .   x2n 


a 2]      f  K    11  n   la  "  '^2    ;     K  n   ,*2 ,   ,Qnv 


(a  -  Yxp  x2 


.Op   is  the  reactor  volume  of  the  2n  th  stage. 
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of  the  system  to  determine  the  number  c: 

les  or  degrees  of  freedom  is  often  re- 
designing and  optimizing  comples  systems.   In 

lysis,  an  N  stage  system  is  considered;  however, 
each  stage  is  assumed  to  consist  of  a  mixing  stage  and  a 
Lon  stage;  thus,  the  actual  model  we  consider  is  ccr 
posed  of  21."  stages.   (see  Figures  7  and  6) 

rid  Number  of  System  Variables 
a)   Volumetric  flow  rate  of  feed, 


ef1"1,     n  -  1,1  ...  N;    N 


b)  Volumetric  flow  rate  through  each  stage, 
xj,   i  =  0,  1,  .. .  2N;   2N  +  1 

c)  Organic  concentrations, 

x*    i  -  0,  1,  ...  2     2N  +  1 

::.   concentrations, 


-  0,  1,  ...  2N;    2N  +  1 

4 
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e)    ..actor  volume  at  each  staj 


r.  =  i,  ...  I.;    N 


)   Concentrations  of  crganics  and  organisms  in  feed, 


xi  and  x; ;    2 


/.  Total  number  of  system  variables  =  3(2N  +  1) 

+  2N  +  2  =  cX  +  3. 

(2)   Type  and  Number  of  Relations  Among  System  Variables, 

a)   Material  balance  of  organic  component, 
(i)   Mixing  stages, 

2n-2   2n-2   _  2n-l  f    2n-i   2n-l 
x-±  x2    +»2    x2  =  x-|_    x2    ,   n  =  1,  ...  N;   N 

(ii)  Reaction  stages  (Tank  reactor  is  under  consider- 
ation) , 


2n-l  2n-l       „  2n    2n  2n        -      M   M 
x-     Xp    +  rg  6  2  ~  x~\      x2  '    XI  =  1,  ...  N;   w 


b)   Material  balance  of  organism  component, 
(i)   Mixing  stages, 


xf-2  xf-2  +  ef"1  x£  =  x2""1  xf"1  ,  n  -  1,  ....  Mj   N 


(ii)  Reaction  stages, 
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-1   v,  2n-l    2n>       .       .. 
=  .       -  x2  J ,   n  -  1 ,  . . . ,  N;   a 

c)   Overall    terial  balance  at  eacr.  _tage 
(i)  s , 


2n-2    ;  .  -1    2n-l       ,       AT 
=  Xn    ,   n  =  1,  . . . ,  L; 


(ii)  Reaction  stages, 


2n-2 
xr-  -  -  xr-  -,   n-1 N;   N 


Total  number  of  relations  =  6N. 

(3)   Degree  of  Freedom  of  System. 

From  the  total  number  of  variables  and  relations  ob- 
ta -     in  previous  sections,  the  number  of  degrees  of 
fre   .  for  the  entire  system  is 

.  -   N  +  5-6N=2N+5. 

e  system  is  specified  by  the  variables,  x.:,  XZ> 
-:'"',  :•:'",  x  ,  xf  ,  then  the  number  of  degrees  of  freedo: 

4    1     *-     ^    4 

F  *  2N  +  5  -  7  -  2N  -  2  . 


->  > 


2.      lation  ^::;dy 

In  cr^cr   to  obtain  a  better  understanding  of  the  b  - 
havior  of  the  system  represented  by  the  proposed  model,  a 
simulation  study  has  been  carried  out.   Because  it  is  dif- 
ficult to  represent  graphically  a  system  with  a  high  number 
of     isions,  only  a  two-stage  completely  mixed  flow  re- 
actor system  has  been  considered. 

The  N  stage  system  is  schematically  shown  in  Figures  6 
and  7.    ..en  N  =  2,  this  diagram  corresponds  to  the  two- 
stage  completely  mixed  flov;  system  used  in  this  simulation 
study.   In  accordance  with  the  derivation  in  the  previous 
section,  two  degrees  of  freedom  exist  in  this  problem 
(i.e.  ?  -  2*2  -  2  =  2).   The  amount  of  the  allocation  of 

wasoe  water  in  the  first  stage,  Or;,  and  the  reactor  volume 

2 

cC   the  first  stage,  Q 2,    are  chosen  as  being  the  two  decision 

variables  for  this  study. 

By  using  the  material  balance  equations,  Equation  (2.8) 
through  (2.20),  as  derived  in  the  previous  chapter,  and  the 
assigned  boundary  conditions,  x-,  ,  x;:,  XT',  xt~,  xt ,  x*  ,  x| , 

■L    <c    if    J.    <~     c.  *f 

and  the  two  specified  decision  variables,  0  ,  G~-,  the  overall 

1    2 

accumulated  reactor  volume,  xtJ,  can  be  evaluated.   The  sim- 
ulation study  has  been  carried  out  by  systematically  altering 

the  values  of  67  and  Q^   and  computing  the  value  of  xo  for 

x  2  -"> 

each  case. 


a  are   employed:      Y  «=  0.5,   k  s     =0.1 

br"1,  K  ■  1       ...,;,--  .       ..     ■        mg/e  ,   x:  -  1800 

hr,    x£  -  1!  4  ,    xj  000  rng/z  ,   x^  =  6000  V'r.r ,    x£  = 

■ 

-  cal  computation  for  this  simulation  has  been 

an  IBM  computer  1620  and  the  numerical  re- 
sults are  listed  in  Table  1. 

re  11  shows  a  plot  of  the  allocation  of  sewage  in 
first  stage  versus  the  overall  accumulated  volume  with 

le  of  the  first  stage  as  parameter.   From  this  figure 
a  two-dimensional  surface  level  diagram  is  in  turn  plotted 
in  Figure  12. 

From  Figures  11  and  12  it  can  be  seen  that  the  objec- 
tive function  for  minimizing  the  overall  reactor  volume  is 
a  ui  .1.1  function.  The  minimum  point  is  located  at  the 
positio]  re  3900  liters  per  hour  of  sewage  feed  are  in- 
troduced into  the  first  stage  and  where  the  reactor  volume 
the  first  stage  is  6250  liters. 

Figure  12  are  also  plotted  two  searching  loci  ob- 
.  discrete  version  of  the  maximum  principle 
teepest  descent  search  technique  as  discussed  in 
hapter.    ie  one  stage  system  can  be  obtained  as  a 
case  in  this  simulation  study  by  assigning  no  reactor 
llocation  of  sewage  feed  in  the  first  stage. 
jcial  case  of  one  completely  mised  reactor  system, 
ctor  volume  of  14,833  liters  is  required  to  obtain  the 


Tabic  1.   Overall  Reactor  Volume  Computed 
Simulation  of  a  System  Compos 
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f   the   first  reactor  gs  parameter. 


30. 


- 


70C 


CO 

V- 

CL> 

.tl 

3 

i_ 

6500 

o 

■«- 

o 

o 

L_ 

6000 

CO 

!— 

«J^ 

CD 

•4- 

5500 

v*— 

o 

G3 

£ 

5000 

J3 

O 

> 

r» 

4500 

_> 

3300         4000 


<-ic-j 


4400 


0,  ,   Allocation     of    sewage    to    the    first   stage  (l/hr.) 


Fig.  12.   Surface    contours   for  the   system  and   the 
searching    locuses. 
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.  .        Accc-         he 
..s  section,  there  is  no  optimization 
...    stem  sir.ce  no  ce~ree  of  freedom  exists. 


>  J 


C.       [V 

o:^  ^s 

In  the  previous  chapter,  the  simulation  study  ha    .en 

^d  for  the  two-stage  complete  mixing  system.    hen  a 

re  complicated  system  is  being  treated,  it  becomes  very 

time  consuming  to  use  simulation  to  optimize  the  system. 

^'r.ere   are  a  number  of  techniques  that  can  be  used  to  op- 

.  Lze  complicated  systems.   Some  of  these  include  (A) 

fferentiation,  (3)  linear  and  nonlinear  programming, 

dynamic  programming,  (D)  experimental  search  techniques, 
.  i  (E)  the  maximum  principle.   The  descrete  maximum  prin- 
ciple is  a  powerful  technique  for  treating  multistage  se- 
mtial  systems,  and  it  will  be  used  to  optimize  several 
models  of  the  system  we  are  considering.   However,  before 
treating  the  multistage  step  aeration  systems  a  brief 

cription  of  the  discrete  maximum  principle  will  be  given. 

rHE  DISCRETE  MAXIMUM  PRINCIPLE  (10) 

(1)   Statement  of  The  Algorithm  for  Simple  Processes. 

A  schematical  representation  of  a  simple  sequential 
Lse  process  is  shown  in  Figure  13.   The  state  of  the 
process  stream  denoted  by  an  s-dimensional  vector,  x  - 

:  ,  x9,  ...,  x  ),  is  transformed  at  each  stage  accordin 
to  the  decisions,  denoted  by  a  t-dimensional  vector,  d= 
,_,    ^j  •••j  ®t)i  r'a^e  at  that  stage.   The  transformation 
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of  the  process  stream  at  the  nth  stage  is  described  by  a 
set  of       lance  equations  represented  cy 

n         n    ,   n-1       n-1  n-1     Qn     _n  n, 

xi  "  Ti    *xl      '    x2      >    '  *  *  »   xs      )   °i)   ^2'    •••0t)>    x  ~  x> 

2,     .  .  .    s  (1) 

•  is  a  transformation  operator  for  component  i  at  the 
nth  stage.   The  optimization  problem  which  we  want  to  con- 
sider here  is  one  in  which  all  the  initial  conditions,  x-,, 
Xp,  ...,  x^j    are  given,  and  one  of  the  end  conditions, 
x",  is  preassigned.   For  this  optimization  problem,  we 

K 

assume  the  objective  function  is  of  the  form, 


S  =  ZI  c,x^  .  (2) 

i=l  x  x 

i^k 


rhe  problem  we  wish  to  consider  is  one  in  which  we  want  to 
extremize  this  objective  function  subject  to  these  given 
initial  and  end  conditions,  and  the  equality  conditions 
/en  by  Equation  (1). 

al  ;orithm  for  solving  such  an  optimization  problem 
the  discrete  maximum  principle  is  to  introduce  an  s- 
mensional  adjoint  vector,  zn,  and  a  Hamiltonian  function 
satisfyin  i 

-  £     zT:   T?  (x11"1;  Sn),   n  =  1,  2,  ...,  N        (3) 
i=l  x  x 


n_] 


'i     .  r.-l 


i  *  1,  2,  ...  s  U) 

n  =  1 ,  2 ,  .  .  .  T." . 


an; 


=  c  j        l  «  1 .  2 ,  . . .  s  1 5  J 

i  f  k  . 


.1  seouence  of  decisions  "   is  determined  from  the 

i 


L  t  i  o  n , 


=  0  ;      1=1,2, ...t  (6) 


n  =  1 ,  2  ,  .  .  .  N  , 

. .  the  decision  variables  lie  in  the  interior  of  the  domain 
Lon  variables.   If  the  optimal  decision  variables 
lie  on  the  boundary  of  the  domain,  then  the  Hamiltonian 

-ion ,  H  ,  should  be  an  extremum  at  the  boundary.   Both 
x  ano  z  are  considered  as  fixed  in  extremizing  the  Hamilton- 
. 

(2)   Computational  Scheme. 

- jai  and  numerical  solutions  will  be  outlined 

.  occ:.     .  Lch  one  end  stage  state  variable  3t  is 

K 


i  -  1,    2, 

•  •  •  «    o 

n  -  1,    2, 

.  .  .  ,     ..     , 

i  -  1,    2, 

•    •    •    j         O 

n  =  I,    2, 

...,    N 

cal  solution.   Referring  to  Figure  13,  a 
litions  x°  are  given,  and  x  is  also  fi: 
„re  i  =  1,  ...,  s).   Now  we  are  to  find  the  :.    Leal 
values  for  the  following  variables  at  the  optimal  con 
tioj 


n 

4  > 


"i 


N 

LACCJO   A.     . 

k 

n  =  1 ,  2 ,  .  .  . ,  N  , 

s 
-  e  total  number  of  unknowns  as  established  above  are 

t«N  +  2-s-K  +  1  -  1  =  t-M  -r  2-s-N. 

The  available  relations  are  summarized  as  follows: 

Relation  Nur.be r  of  relations 

1.  Performance  relations,  Equation  (1)      s'N 

2.  =  Cj_,  i  f  k,  Equation  (5)  s-1 

3.  Adjoint  relations,  Equation  (4)        s'(N-l) 

•:n 
A-*   :lj—  =  0,  Equation  (6)  f  N 

5.   Djective  function  Equation  (2)  1 


u 


refore, 
s'N  +  s'(N-l)  +  (s-1)  +  (t'N)  +  1  =  f  a  t   2-:  •:;. 

.  r.ce  the  .    er  of  variables  is  equal  to  the  numt  .r   of 
relations,  the  desired  quantities  can  be  found  by  solving 
e  ecu.      .     .  ar.eously.   The  simultaneous  solution  cf 
,se  equations  is  extremely  difficult  except  for  some 
Le  cases.    .cause  of  this,  numerical  procedures  fcr 
treating  the  particular  type  of  problems  under  consideration 
re   been  developed  (10). 

srical  solution 

.  ....  rical  computation  for  the  process  can  be  car- 
ried  out  by  the  following  steps. 

Step  -.   Assume  a  set  of  values  for  9?  (i  =  i,  2,  ..., 

N  N  N 

t;  n  =  1,  2,  ...,  N)  except  for  6^.   Since  x,"  is  given,  d. 

11  be  a  dependent  variable. 

2.   Calculate  all  the  state  variables  x?(n  =  1, 
2,  ...,  N-l;  i  ■  1,  2,  ...,  s)  by  using  the  correspond!. 
mce  equations. 
tp  3 «   Calculate  Q  ■"   from  the  performance  eouaticn, 


'i 

frc  performance  equation,  i.e., 


.  .   Calculate  x     B  1,  2,  ...,  k-1,  k+1,  ...,  s) 


x\   =  T^    (x        ,  6  )  ;  i  =  1,  2,  ...,  k-1,  k+1,  ...,  s, 

value  obtained  from  steo  3- 
k 


so  5.  Calculate 


dx. 

— - —  and  — - 


with  all  state  variables  and  decision  variables  obtained 
from  previous  steps. 

In  a  multistage  operation  with  one  end  stage  state 

variable,  say  xjj,  fixed,  the  adjoint  variable  associated 

N 
with  it,  namely  zj.  should  be  dependently  determined  to 

hold 

auN   i.   .rax? 

^r-=  Z     z)   -|  =  0  ,  (8) 

"k     1  =  1      S9k 

thus  obtain 


s     dxN 

-Z  zN   i 
i-l  1  _ 

=  i£ ^   .  (9) 

d: 

k 


Step  6.  U^    Equation  (5)  and  the  corresponding  I  - 
rivatives  obtainea  from  Step  5,  zg  can  be  calculated  from 
iion  (9) . 


_ 


Calculate  the  r    .  .     Dint  \       s, 
■1,  2,  ...,  N  -  1;  i  =  1,  2,  ...,s)  from 

•  .  ) . 

Calculate  all  values  of  ■'■      (n  =  1,  2,  ..., 

dG"- 

i  =  1,  2,  ...,  t)  except froi        n  (6). 

•  -:: 

i  ues  of  t>  are  not  the  optimal  values,  the 

3  O 


f  ]  (n  =  1,  2,  ...,  N;  i  =  1,  2,  ...,  t)  will  not 


be  zero.   We  define 


iep  9.  Using  the  method  of  steepest  descend  a  new 
set  of  decisions  can  be  found  to  replace  the  old  set  ac- 
cording  to  the  following  equations 

<G'.    .  =  (ei\ld  +^Gi>    i  -  1,  2,  ...,  t        (10) 

n  =  1,  2,  .  .  .  ,  N 

i  =  1,  2,  ...,  t;  n  =  1,  2,  ...,  N  are  specified 
_  corr         _ues  whose  sign  is  determined  by  that  of 

•  •  !  ^i« 
.   If  all  the  values  of  <f>i(n  -  1,  2,  ... 

...  t)  are  smaller  than  some  preassigned  small 

j  ,  t]       set  of  9  ?  is  approximately  the  optimal 

ration  shou    ..  stopped.  If  not, , 

:„c.-:  to  steo  1. 
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2.          ATI ON  STUDIES  OF  STEP  AERAT] 
OF  SI L  COr     LY-STIRRED  R] 


As  mentioned  in  Chapter  1,  the  distributer.     tie  ar- 
z   load  cr.o  the  flow  behavior  are  i .'.-.portent,  variabl 
in  the  step  aeration  process  which  affect  the  optimal  cc 

of  a  orocess.   In  this  section  the  effect  of  these  variables 

or.       a  required  for  a  specified  degree  of  treatment 

will  be  considered  for  the  step  aeration  process  composed 
of  several  completely-stirred  tanks  in  series. 

he  problem  is  to  determine  the  distribution  of  the 
organic  load  which  minimizes  the  total  required  reactor 
volume  for  such  a  system.   In  Chapter  2,  this  system  is 
described  in  detail  and  the  related  system  equations  are 
rived  (see  Figures  6,  7,  and  8). 

(1)   Formulation  of  the  Problem 

problem  is  oo  find  the  optimal  allocation  of  the 
sewage  and  the  optimal  distribution  of  volumes  between  al- 
locations for  carrying  out  the  waste  treatment  process  in 
a  sequence  of  completely  stirred  tank  reactors  so  that  the 
total  reactor  volume  is  minimized.   The  objective  function  o: 
this  problem  is  therefore  simply  represented  as 

S  -  x  (II) 

j 

2  v 
.re  xV"  represents  the  total  reactor  volume. 
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latior     re   be^. 

shown  -         6,  7  and  8  the  od 

E  account  -     process  which  is  described 

as  follows: 

-  ^uation  {2.6),   the  overall  material  balance  is 

-1  =e-":  -r   x?n-2  .  (12) 

1-1 

?rc     Lation  (2.10),  the  material  balance  of  the  oxidizable 
nics  is 

a2n-l       2n-2  2n-2 

,--  -  L ^  T  ^  .  x-  (13) 

Q2n-1  ^n-<i 

ation    (2.13),    the  accumulated   volume   is 


•  -1  ---2  M  .  » 

=  x^  u4) 


e  ever,  numbered  stages  account  for  the  biological 

ich  is  described  mathematically  as  folic  : 

.  (2.14),  the  overall  material  balance  leads  to 


(15) 


ns  (2.20)  and  (2.21),  the  material  I    nee  of 
nics  gives 


- 


s    2n  [  2n-l         *2 
2n-l   .--    2n-l       A2n  "   ;  z/  ^  x: 


^ 


;,  -  -2 


=  0  .  (16) 

i  ..   16)  can  be  generally  expressed  as 

2n   p  /  2n-l  v2n-l  ~2n^  r7l 

=  i  (x^    ,  X^    ,  6  ^  J .  (17) 

From  Equation  (2.19),  the  accumulated  volume  is 


2n-l   a  2n  /-  ■  * 

3     3       2   ' 


In  applying  a  discrete  form  of  the  maximum  principle  to 
determine  the  optimal  policy  for  the  system  as  stated 
above,  the  following  Hamiltonian  functions  and  adjoint 
variables  are  introduced.   The  Hamiltonian  functions  for 
the  cing  stages  and  reaction  stages,  respectively,  are 


H2n-1  =  z  z2n-l  x2n-l^   n  =  1,  2,  ...,  N,  (19) 

i=l  x     x 

H2n  =  £     z?n  x?n,   n  =  1,  2,  ...,  N  .  (20) 

-1  x        x 


i   adjoint  variables  for  the  mixing  and  reaction  stages, 
respectively,  are 


J 


.  - —  / 


I,  .  -  .    .0 

-     Lj    ...    s  -.^    3 j    respe        /ely 
.-1 


"•"     - 

-1 

>        J 

• 

! 

,_ 

>  «-^  . 


[24] 


—  - -    .  - -^-  j 


5 


a  2n  v2r.     r.ax  o  |  fo 

z2n-l  .  z2n     ^1    (K  +  xf )    ,  (27) 


2n-l  2n 

.3  -3      , 


(28) 


where 


_,       v   2n-l        0   2r.-l     2n       ^  2n  ,  s  2n  2n-l     2n-l 

B  =  Kx1  +  2x1  *2      ~  *  2     ^  x2      -  xx  x£ 


fl2n    wrr.ax    /      "o      ,    Y     f       v     2nx 
D2     T^   (~2l~i  +  Y  x2  -  Y  x2    K 

xl 


Since  x^*'  and  x2  are  specified,  the  only  boundary  condition 
for  the  adjoint  variables  in  this  system  is 

23'"  -  1.  (29) 


From  Equations  (25),  (25),  and  (29),  we  thus  obtain 


=  z^n  =  1,    n  =  1,  2,  ...,  ...  C 


so-Cc ik   for)  he   discr* 

optical   decision   va.  ich 

stc  the    equatic. 

^     •  n-i  2 =0  n  =  2  K  (31) 

l— _  c  t 


3      0    ax2n 

-  =    Z      zfn  — f—  =  0    ,      n  =  1,    2,    ...,    N.  (32) 

1-1  ,c2 


.  1  decision  variables  are  in  the  interior  of 
differentiating  Equations  (12),  (13) >  and 
(14)  ■   .  respect  to  61  "  ,  Equation  (3D  becomes 


0    _ 

-    2n-l   r2r-l 

2n-l  "  2i   —2 

(ef-1  ♦  .   •"-2)x;-(6f-14  +  xf-2  xf-1) 

n  =  1,  2,  ...,  N  .  (33) 

rig  Equations  (15),  (16),  and  (18)  with  re- 
el tc       i  substituting  the  resulting  expressions  into 

„  orr.es 


>> 


2n      0  x%n   k°     A         r 

z2n  _2 x  (__o v  f  _  v   .       .  0 

2H   ~  22    B   Y    lY2n-l    *  X2    r  X2  ;    -1   U  ' 

1 

-  i,  2,  ... ,  ...  (: 

If  the  optimal  decision  variables  are  on  the  boundary 
the  domain,  the  corresponding  value  of  H  1_  cr  H  "  must 
be  an  extremum. 

(2)   Computational  Scheme 

. ..  rical  solutions  to  the  system  of  algebraic  equation* 
given  by  Equations  (12),  (13),  (14),  (15),  (16),  (IS),  (23), 
(24),  (26),  (27),  (30),  (33),  and  (34)  are  obtained  by  fol- 
lowing the  computational  scheme  presented  previously.   For 
this  problem,  the  step  by  step  computational  procedure  is 
as  follows: 

Step  1.   Assume  values  for  6?n~  ,  e?n,  n  -   1,  2,  ..., 

..-I. 

0   0   o 
Step  2.   Using  the  given  conditions,  i.e.,  x-,,  x~ ,  x~  , 

x.  ,  x   ,  calculate  the  state  variables,  x~_    ,  ant  xr  , 

i  =  1,  2,  3,  n  =  1,  2,  ...,  N,  and  ~   ,   «  from  Equations 

(12;  through  (IS). 

3p  3-   Using  Equation  (34)  for  n  =  N,  calculate  z£  . 

_  - 
Step  4.   Using  Equation  (27)  for  n  =  N,  calculate  z*   . 


54 

for  n  =  lv,    calcu]  -1 

......   uj  ...      -  s    (23,     (20,     (26),    ^r.e    (27), 

rest  of  the  ac joint  variables  z^n-l^  zfn~-, 


;n,  n  =  1,  2,  .  .  .  ,  N-l. 


-1      ,2n 


p  7.   Calculate  the  values  of  j;-/    ar 


o62 


-5  (33)  and  (34),  which  are  represented  by  6     ana 

On  "1         On 

.  !      .ely.   (Note:   If  the  values  of  t*:*'~x   and  Q~ 

_2 

art.      ie  optimal  values,  the  values  of  fytr-~~   and   n  will 

^  zero) . 

Step     Using  the  following  equations,  calculate  the 
.  set  of  decision  variables. 


[n_1'new-  t9in"1)old+  (-e^Xsgn^f-1) 


<61    ..=   ^f'olc-  t^Hsgn^f) 


_,  2 ,  ...,  ..  —  _ 


-  1  if 

-1      ^  0 


-/> 


iQy1   are  specif iea  small  constants. 

St  en  9.     all  the  values  of  ^i^"1  and  ofr-    (n  =  1, 

. ..,  .-__.■  are  .     sr  than  some  preassi  ..   small  va 
.,      the  set  of  6  jr11"1  and  G^n  (n  =  1,  2,  ...,  N-l)  i 

~ely  the  optimal  decision  set  and  the  iteration  ca 
stopped.   If  not,  the  computation  should  go  back  to  Step 
1. 

In  the  present  work  the  numerical  computation  according 
to  the  above  procedure  has  been  performed  with  an  IBM  1620 
computer.   The  computer  source  program  for  a  two-stage 
system  is  given  in  Appendix  1. 

A  .  lj  h  modification  of  the  above  numerical  procedure 
is  required  for  Lhose  systems  in  which  there  is  no  allocation 
to  one  cr  more  stages.   The  inequality  constraint 


e^n-i   >y  0 


be  satisfied  at  each  stage;  that  is,  the  allocation  to 
each  tar./:  must  be  either  positive  or  zero.   This  inequality 
istraint  is  treated  by  setting  d?n~x   equal  to  zero  for 
;e  in  which  a  negative  value  for  9?""1  is  cbtai: 

-cal  computation.   For     2  stages  with  zero 
allocation,  9  2n-l  should  be  positive  ana  H~r-~-  should  be  a 
at  the  optimum.   Thus,  the  zero  value  for  &  f n~"L  is 
■  in  the  ..   rical  computation  as  long  a.     ~~  is 
p  o  s    e . 


.  6  .  cussi 

•suits  '.he  optimal  i. 

Lmal  vo]         bution  for  the 
.  r.  system  consistin  ;  of  se    -  completely 

s  in  series  are  presented  in  this  section, 
tic  constants  are  assigned  the  foil      valu 
,  K  =  100  mg/liter,  ieVv  =  0.1  hr"1.   The 

follov ditional  conditions  are  also  specified:   ..  - 

0.4,  :■-   -  150  mg/4  ,  x^  -  8000  nsg/j  ,  x^  -  800  m&k   ,  x*1  = 


"1=  4500  tyhr. 


rical  solutions  have  been  obtained  using  these 
for  steo  aeration  systems  containing  one,  two,  three, 
four,  and  five  completely  stirred  tanks-in-series.   The 
.  .    i  volume,  the  optimal  influent  flew  rates,  ar.o  the 
volumes  are  presented  in  Table  2  for  each 
le  relationship  between  the  minimum  volume  of 
3  sys  .  mber  of  completely  stirred  tanks 

system  is  shown  in  Figure  14.   he  curve 

Le  set  of  discrete  points  shows  that 
red  aeration  system  volume  decreases'  exponen- 
er  of  tanks  in  the  system  is  increased. 
:.t  effects,  the  extent  of  longi- 
.  eration  sy.      d  the  method  of 
to  the  system. 


...... 


^  13000 

3  12000 

z 

§  nooo 


o 

5- 


ICOOO 


^iil 


__- 


12  4  5  6  7 

The  number   of  completely  stirred    tanks. 

Fig.  14.   Optimal   total  reactor    volume    v.s.  the 

number  of  operating  tanks  for  completely 
stirred  tanks   connected   in  series. 


Lon  sy 

pletely  £      I  tc s  con- 

e  s  in  j  one  stage  system  repre- 

s  a  system  with  comp] 

ree,  four,  and  five  -     systems  correspond  to 
systems  with  less  longitudinal  mi    ;.   As 
.  of  cc.       -  stirre        in  the  systerr.  increases 
mount  of  1;      linal  :..     ;  decreases.   Figure  14 

t  the  required  volume  for  the  aeration  systerr.  de- 
creases  as  the  number  of  tanks  in  the  system  increases.   This 

hat  as  th^     nt  of  longitudinal  mixing  is  ae- 
creased,  the  required  volume  decreases.   The  results  for  the 
ree,  four,  and  five  stage  systems  in  Table  2  show  con- 
at  Ion  ;itudinal  mixing  is  undesirable  in  the 
■art  of  the  aeration  system  (the  optimal  volume  of 
ould  be  zero  if  longitudinal  mix!    -re 
desirable).         suit  agrees  with  that  of  Bischoff  (11) 
own  that  longitudinal  mixing  is  often  undesirable 
e  last  part  of  an  aeration  system. 

results  in  Table  2  also  show  that  the  optimal  mcth 
influent  to  the  aeration  system  involves  al- 

ore  than  one  location.   The  Michaelis* 

.-ation, 
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is  used  in  this  work  to  indicate  that  the  growth  rate  de- 
r.as  or.  roth  the  concentration  of  organics  and  organisms. 
an  see  two  counteracting  effects  that  influence  the 
growth  rate  when  influent  is  added.   The  addition  of  in- 

Luent  increases  the  concentration  of  the  organic  substrate 
and  decreases  the  concentration  of  the  microorganisms.   As 
shown  in  Table  2  the  optimal  method  of  feeding  the  five 
stage  step  aeration  system  involves  distributing  most  of 
the  influent  to  the  first  half  of  the  system.   In  the  first 
part  of  the  system  it  is  desirable  to  maintain  as  high  a 
growth  rate  as  possible  by  addint  the  influent  in  an  op- 
timal manner. 

Influent  should  not  be  fed  to  the  last  part  of  the 
aeration  system,  because  it  is  necessary  to  reduce  the  or- 
ganic concentration  to  a  specified  small  value  at  that  end 
of  the  system.   Near  the  exit,  the  addition  of  influent 
would  increase  the  growth  rate,  but  it  would  also  greatly 
increase  the  organic  concentration  which  has  to  be  reduced 
to  a  specified  value.   As  shown  in  Table  2,  the  optimal  re- 
sult is  not,  to  feed  any  influent  to  the  last  part  of  the 
aeration  system. 

The  effect  of  the  organism  concentration,  the  organic 
icentration  in  the  return  sludge,  and  the  return  sludge 
w  rate  on  the  optimal  design  have  been  studied  for  the 
ration  system  containing  two  completely  stirred  tanks  in 
series. 
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In  Figure  15  and  Table  3  are  shown  the  optimal  policies 
'  the  system  for  various  inlet  organism  concentrations.   The 

curve  which  relates  the  total  optimal  volume  to  the  inlet 
organism  concentration  shows  that  an  increase  in  the  organ- 
ism concentration  in  the  return  sludge  will  reduce  the  total 
required  volume.   This  is  due  to  the  effect  that  an  increase 
in  organism  concentration  will  favor  the  rate  of  growth. 
Curves  which  relate  the  optimal  influent  distribution  and 
the  optimal  volume  distribution  to  the  inlet  organism  con- 
centration are  also  shown  in  Figure  15.   These  results  show 
that  as  the  inlet  organism  concentration  is  increased  a 
greater  portion  of  the  influent  to  the  system  should  be  fed 
to  the  second  stage  and  a  greater  portion  of  the  total  volume 
should  be  allocated  to  that' stage.   Except  for  the  organism 
concentration  the  data  used  in  this  study  are  the  same  as 
those  given  at  the  beginning  of  this  section. 

Similarly,  the  optimal  policies  for  the  two  tank  step 
aeration  system  with  various  inlet  organic  concentrations 
in  the  return  sludge  have  been  studied  as  shown  in  Figure  16  . 
and  Table  4.   However,  the  data  used  in  this  part  of  the 
present  study  are  slightly  different,  viz.,  Y  =  0.5  mg/mg, 

K  =  100  mg/l.,  kgax  =  0.13  hr"1,  R  =  0.334,  x£  =  4000  mg/l., 
x*  =  800  mg/l.   Figure  16  shows  that  the  total  optimal  volume 

is  nearly  proportional  to  the  organic  concentration  in  the 
return  sludge.   The  curves  which  relate  the  optimal  influent 
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ution  and  the  optical  volume  distribution  to  the  inlet 
or..":  Lc  concentration  indicate  that  as  the  inlet  organic 
concentration  is  increased  a  greater  portion  of  the  influent 
zc   the  system  should  be  fed  to  the  second  stage  of  the 
systeir.  and  a  greater  portion  cf  the  total  volume  should  be 
allocated  to  that  stage. 

These  results  imply  that  an  increase  in  either  the 
~anism  or  organic  concentration  in  the  return  sludge  will 
result  in  a  smaller  fraction  of  the  influent  and  volume  be- 
ing allocated  to  the  first  stage.   This  effect  appears  to 
be  due  to  an  effort  to  take  advantage  of  this  increased 
inlet  concentration  in  the  first  stage  by  reducing  the 
amount  of  influent  added  to  that  stage. 

Increasing  the  recycle  ratio  will  reduce  the  volume 
required  for  the  aeration  system  and  thus  reduce  the  capital 
cost  associated  with  the  volume  requirements  of  the  system. 
Hcwever,  the  costs  associated  with  the  sedimentation  process 
ar.d  the  cost  of  pumping  the  recycle  stream  will  increase 
with  an  increase  in  the  recycle  flow  rate.   Furthermore,  as 
ir.cre   sludge  is  recycled  through  the  system,  the  viscosity 

the  aeration  system  will  increase  and  the  cost  of  aeration 
will  increase  because  of  this.   In  order  to  show  that  there 
is  an  optimum  recycle  ratio,  we  assume  a  very  simple  cost 
mcael  in  which  the  total  cost  is  given  by 


CT  -  c1V  +  Z.500  C2R 
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where 

C^  =  unit  cost  of  volume,  $/liter 

hr  / 
C2  =  unit  cost  of  recycle  ratio,  v   /liter 

V  -  total  volume,  liters 

R  ■»  recycle  ratio,  dimensionless . 

In  this  expression  each  value  of  the  total  volume,  V, 
is  the  optimal  value  determined  by  use  of  the  maximum 
principle  for  a  certain  value  of  the  recycle  ratio,  R. 
The  optimal  recycle  ratio  can  be  found  graphically  by 
plotting  the  total  cost  against  the  parameter,  R. 

The  two  stage  model  of  the  aeration  system  has  been 
studied  using  this  cost  model  for  the  case  in  which 
(^/(^  =  0.5. 

In  this  study,  the  following  data  have  been  used. 
Y  =  0.5  mg/mg,  K  -  100  mg/l,  kf>   =  0.1  hr"1,  x§  =  150 

mg/l,  x£  =  8000  mg/l  ,  x£  =  800  mgA  ,  x^N  =  80  mg/l, 

;; 

2£-@2n-l  =  4500  1/hr.   The  results  which  are  presented  in 
n-1  1 

Table  5  and  Figure  17  show  that  there  is  a  minimum  cost  at 

R  ■=  0.55. 

Similarly,  the  optimal  feed  and  volume  distribution 

ratios  are  plotted  in  Figure  18.   The  curves  show  that  the 

ratios  decrease  as  the  recycle  ratio  increases. 
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3.    OPTIMIZATION  STUDIES  OF  STEP  AERATION  SYSTEMS  CONSISTING 
OF  SEVERAL  PLUG  FLOW  REACTORS  IN  SERIES 


As  mentioned  previously,  the  method  of  feeding  the  in- 
fluent to  the  aeration  system  and  the  flow  behavior  are  im- 
portant variables  in  the  step  aeration  process  which  affect 
the  optimal  design  of  the  process.   In  this  section  an 
aeration  system  with  no  longitudinal  mixing  is  considered 
and  the  optimal  distribution  of  the  organic  load  along  the 
aeration  system  is  examined  with  respect  to  its  effect  on 
the  required  aeration  system  volume.   Since  the  undesir- 
ability  of  longitudinal  mixing  in  the  latter  portion  of  the 
aeration  system  has  already  been  pointed  out,  it  is  desir- 
able to  consider  aeration  systems  with  no  longitudinal 
mixing. 

As  in  the  previous  study,  the  system  will  be  divided 
into  mixing  stages  and  reaction  stages  shown  in  Figures 
6,  7  and  8  of  Chapter  2;  however,  this  time  plug  flow  is 
assumed  in  the  reaction  stages.   The  material  balance 
equations  which  describe  the  mixing  and  the  growth  process 
have  been  presented  in  Chapter  2. 

(1)   Formulation  of  the  Problem 

The  problem  is  to  find  the  optimal  allocation  of  in- 
fluent and  the  optimal  distribution  of  volume  between  al- 
locations for  carrying  out  the  waste  treatment  process  in 
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aeration  system  represented  by  a  sequence  of  plug  flow 
reactors  so  that  the  total  aeration  system  volume  is  mini- 
mized.  The  objective  function  of  this  problem  can  be  rep- 
resented as 

S  =  x,  (35) 

where  x.   represents  the  total  aeration  system  volume. 

The  mathematical  equations  which  describe  this  system 
are  almost  the  same  as  those  of  the  continuous  stirred  tank 
reactor  in  series  system  except  for  the  reaction  sections. 
The  performance  equations  for  the  mixing  stages  are  iden- 
tical with  those  mentioned  in  the  previous  sections,  that 
is,  Equations  (12),  (13),  and  (14).   The  performance 
equations  for  the  reaction  stages  are  established  as  fol- 
lows:  From  Equation  (2.15),  the  overall  material  balance 
gives 


x2n  =  x2^1  .  (36) 


From  Equation  (2.30),  the  material  balance  of  the  oxidizable 
organics  gives 


2n  r  k         1 


(a   -   Yx2n_1)  x2n 

X2         ;         K         n     ,X2 


•SB  -  -  l  fe +  fe  «» : — 7te—  -  to  *»<-fci>        '37) 


(a   -   Yx2    )  Xp 
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Equation  (37)  can  be  functionally  expressed  as 

2n    -  ,  2n-l   2n-l  fl2n,  /.^x 

From  Equation  (IS),  the  accumulated  volume  is 

2n    2n-l  &  2n  mqi 

x^   =  x^      ^2   ' 

and  thus,  the  objective  function  is 

.2N 


9   =   Y 
O  —  A 


3 


(40) 


The  optimization  problem  is  thus  one  in  which  we  wish 
to  choose  the  allocation  to  each  mixing  stage  and  the  volume 
of  each  reaction  stage  so  that  we  minimize  Equation  (40) 
subject  to  the  equality  constraints  given  by  Equations  (12), 
(13),  and  (14)  for  the  mixing  stages  and  Equations  (36), 
(37),  and  (39)  for  the  reaction  stages.   It  is  assumed  that 
the  flow  rates  and  concentrations  of  organics  and  organisms 
are  known  for  the  sludge  recycle  and  influent  streams, 
and  that  the  outlet  concentration  of  the  organics  is  fixed. 

The  solution  to  this  optimization  problem  can  also  be 
obtained  by  using  the  discrete  maximum  principle.   As  in  the 
previous  problem,  we  introduce  adjoint  variables  and  a  Kamil- 
tonian  function  for  each  stage.   The  Hamiltonian  function 
and  adjoint  variables  for  mixing  stages  are  given  in  Equations 
(19),  (21),  (23),  (24),  and  (25).   The  Hamiltonian  function 
for  the  reaction  stages  is 
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„2n  m     2n  2n 
zl     xl 


2n  2n 
z2  X2 


z2n  x2n 

z3  x3 


.re  xj*  ,  x~  ,  and  x   are  given  by  Equations  (36),  (38), 
and  (39)  respectively.   For  the  reaction  stages,  the  adjoint 
variables  can  be  obtained  by  diferentiating  Equations  (36), 
(37),  and  (39)  with  respect  to  Xj11"1  for  j  =  1,  2,  3,  and 
substituting  the  results  into  Equation  (22).   Thus  the  ad- 
joint variables  are  given  as 


2r.-l    2r.   1  f   AQ(a  *  kY)  (x^1  -  xf ) 

3  Uclxf-1)2  (a  -  Yx|n)fe-  Yxf"1) 


=  z 


K' 


2n-l 


f   ..2n-l 


(2A   +  Yxt  x, 

o     2   1 


,r.ax 


K^"  (A 


Y*|  xf -1)2 


Jin 


x2n(a  -  Yx2/1"1) 
xf-1(a  -  Yx2n) 


1 


K 


max 


in 


a  -  Yx 


2n-l 


a  - 


Yx^n 


(41) 


where 


K 


a  +  KY 


acx 


2n 


ac(a  -  Yxp  ) 


(42) 


2n-l    2n  1 
2      z2   B 


a  ^  KY 


-T 


K 


/    v  2n  v      2n-l 
ac(a  -  Yxp  )   acx2 


(43) 
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zf"1  -  z\n   .  (44) 

Equations  (25)  and  (44)  reduce  to  z^11"1  =  z^n  -  1,  n  -  1, 
2,  . ..,  N,  since  according  to  the  maximum  principle  the 
boundary  conditon  for  this  adjoint  variable  is  given  as 
z|N  =  i. 

The  necessary  conditions  for  a  stationary  point  can  be 
obtained  by  using  the  so-called  weak  form  of  the  discrete 
maximum  principle  as  given  by  Equation  (6).   For  the  mixing 
stages,  this  gives  Equation  (33)  which  has  been  obtained 
previously;  however  for  the  reaction  stages,  differenti- 
ating Equation  (36),  (37),  and  (39)  with  respect  to  G |n 
and  substituting  the  results  into  Equation  (32)  gives 

?r\  2n 

an2*    z2 


ae|n    B 


+1=0,   n  =  1,  2,  .. . ,  N.  (45) 


(2)   Computational  Scheme 

The  computational  scheme  for  obtaining  numerical  solu- 
tions which  was  discussed  previously  can  also  be  used  to 
solve  this  set  of  nonlinear  algebraic  equations.   However, 
the  step  in  which  we  solve  for  Xpn  is  no  longer  straight- 
forward and  a  numerical  method  such  as  the  Regula  Falsi 
method  or  Bolzano's  method  must  be  used.   In  this  work,  the 
Regula  Falsi  method  was  used.   The  computer  program  for  the 
numerical  solutions  is  given  in  Appendix  2. 
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(3)   Results  ana  Discussion 

Numerical  results  which  give  the  optimal  influent  dis- 
tribution and  the  optimal  volume  distribution  for  a  step 
aeration  system  consisting  of  several  plug  flow  reactors 
connected  in  series  are  presented  in  this  section.   A 
simple  approach  to  attain  a  best  continuous  influent  dis- 
tribution along  a  tublar  reactor  is  qualitatively  considered 
and  discussed  in  the  last  portion  of  this  section.   In  this 
numerical  study,  the  kinetic  constants  are  assigned  the 
following  values:   Y  =  0.5  mg/mg,  K  ■  100  mg/l,  k^ax  = 
0.1  hr~  .   The  following  additional  conditions  are  also 

specified:   R  =  0.4,  x^  =  150  mg/l,  x°  =  8000  mg/l,  x^  = 

SCO  mg/l,  xf*   =  SO  mg/l,  S  G  l*'1   =  4500  l/hr.   Numerical 

n=l 

solutions  have  been  obtained  using  these  data  for  step 
aeration  systems  consisting  of  one,  two  and  three  plug 
reactors  connected  in  series.   The  results  are  listed  in 
Table  6  and  plotted  in  Figure  19.   The  lower  curve  which 
represents  the  results  for  the  plug  flow  systems  shows 
that  the  total  required  aeration  system  volume  decreases 

onentially  as  the  number  of  feed  locations  in  the  system 
is  increased.   However,  the  reduction  in  volume  obtained  by 
adding  other  feed  locations  is  small  for  cases  with  more 
than  two  feed  locations.   This  shows  that  the  effect  of 

influent  distribution  on  the  required  reactor  volume  is 
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small  for  systems  without  longitudinal  mixing.   In  Figure 
19  the  curve  presented  previously  for  the  C.S.T.R.  system 
is  also  plotted  for  comparison  with  the  curve  for  the  plug 
flow  system.   The  two  systems  represent  two  extreme  cases 
in  that  complete  mixing  is  assumed  in  each  tank  of  the 
C.S.T.R.  in  series  system,  while  no  longitudinal  mixing  is 
assumed  for  the  plug  flow  reactor  system.   Other  systems 
with  incomplete  mixing  will  lie  between  these  two  curves 
as  shown  in  Figure  19. 

In  systems  of  plug  flow  reactors  connected  in  series, 
the  volume  of  each  reactor  and  the  quantity  of  influent 
added  at  the  front  of  each  reactor  does  not  influence  the 
longitudinal  mixing  in  any  way.   On  the  other  hand,  in  a 
system  composed  of  a  series  of  completely  mixed  tanks  a 
change  in  the  number  of  tanks  or  even  the  volume  of  one 
tank  will  change  the  mixing  pattern  of  the  system.   In 
Figure  19  the  decrease  of  required  reactor  volume  with 
increasing  stage  number  for  the  plug  flow  system  is  due 
only  to  the  effect  of  changing  the  feeding  pattern  for  the 
influent  entering  the  system,  while  for  the  C.S.T.R  in 
series  system  the  decrease  is  due  both  to  changes  in  the 
mixing  pattern  and  to  changes  in  the  influent  feeding 
pattern.   Another  graphical  representation  of  these  systems 
can  also  be  made,  and  it  is  shown  in  Figure  20.   Tho  solid 
curves  shown  in  Figure  20  are  a  plot  of  the  reciprocal  of 
the  reaction  rate  versus  the  organic  concentration  with  the 
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Fig.  20.    Growth  rate   trajectory   in  an    optimal   three   stage 
plug    flow    step    aeration     system. 
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flow  rate  as  a  parameter.   The  equation  of  the  solid  lines 
is  obtained  from  equation  (2.18),  which  is  as  follows 

(K  +  xj 
I       i .  (46) 

kmax*2   ^  +  Y4-  Yx2> 


where 


A   .00   v00   v  0  f 

Aq    —  X-i  Xi   t  I  X-i  X^>  ™  I  X-i  Xp  . 


The  dotted  lines  are  determined  from  the  mixing  relation, 
equation  (2.10),  which  can  be  written  in  the  form 

Y2n-2  /vf   Y2n-2} 

x2n-l  =  xf  .  i U2  .  X2 I  (4?) 

x2      x2        v2n-l  ^M 

xl 


In  Figure  20  the  heavy  solid  line  shows  the  path  of 
the  optimal  locus  of  the  three  stage  plug  flow  reactor  in 
series  system.   The  recycle  stream  is  at  point  A ,  x j_  =  1800 
and  Xp  =  150.   In  the  first  stage  where  1050  liters  of  in- 
fluent are  added  to  the  recycle  stream,  curve  AB  up  to  point 
B  is  followed.   As  this  material  passes  through  the  first 
reaction  stage,  the  reaction  curve  BC  is  followed  to  point 
C.   The  area  enclosed  by  bcCB  represents  the  holding  time 
required  for  the  first  reactor.   In  the  same  manner,  the 
curve  is  followed  up  to  D  as  2500  liters  of  the  second  in- 
fluent are  added  to  the  second  mixing  stage  and  then  up  to 
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j.  21.   Growth  rate    trajectory   in   an   optimal   five   stage 
completely    stirred    tank    step  aeration    system- 
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E  as  the  rr.aterial  passes  through  the  second  reaction  stage. 
Similarly  for  the  last  stage,  the  curve  is  followed  up  to 
F  and  then  to  G.   Similarly,  the  areas  enclosed  by  dDEe  and 
fFGg  represent  the  holding  time  required  for  the  second  and 
third  reactor  respectively.   The  total  holding  time  is  the 
sum  of  these  three  areas. 

By  using  this  graphical  representation,  the  advantage 
of  employing  a  plug  flow  reactor  rather  than  a  completely 
stirred  tank  flow  reactor  can  be  clearly  illustrated.   Based 
on  the  graph  used  in  Figure  20,  the  optimal  locus  of  the  five 
completely  stirred  tanks  in  series  system  is  plotted  in  Fig- 
ure 21,  i.e.  the  curve  ABB'CDD 'EFF ' GG'HH'I .   The  total  hold- 
ing time  will  be  the  sum  of  the  five  rectangular  areas;  they 
are  cbB'C  edD'E,  gfF'G,  hgG'H  and  ihH'I.   For  the  same  in- 
fluent flow  situation  and  a  plug  flow  reactor,  the  locus 
will  be  along  ABCDFGKI.   The  areas,  BCBT ,  EDD • ,  GFF' ,  KGGT , 
and  IHHT,  represent  the  amount  the  holding  time  is  reduced 
by  using  plug  flow  reactors  instead  of  stirred  tank  reactors. 

From  a  graphical  analysis,  we  can  see  that  a  continuous 
feed  allocation  along  a  plug  flow  reactor  which  keeps  the 
organic  concentration  at  200  mg/liter  will  yield  the  opti- 
mal policy  for  this  system.   In  other  words,  the  optimal 
locus  will  be  the  locus  ABCD  passing  through  the  minimum 
points  of  the  mixing  lines  as  shown  in  Figure  22.   A  nu- 
merical result  has  been  obtained  by  using  a  finite  approx- 
imation of  the  continuous  allocation  system.   In  the 
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22.  Growth  rate  trajectory  in  a   plug   flow    step  aeration 
system    with  optimal    allocation   along    the   length 
of    the    system. 
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approximate  system,  influent  is  added  at  each  mixing  stage 
until  an  organic  concentration  of  200  mg/liter  is  reached. 
A  2%   BOD  reduction  in  each  reaction  stage  is  allowed;  that 
is,  the  holding  time  is  chosen  so  that  the  exit  concentra- 
tion from  each  plug  flow  reactor  is  196  mg/liter.   This 
pattern  is  followed  until  all  of  the  influent  has  been  added 
to  the  system.   The  total  required  volume  by  this  approxi- 
mate calculation  is  8025.2  liters.   This  is  an  approximation 
to  the  limiting  value  of  the  two  curves  shown  in  Figure  18. 
When  the  number  of  tanks  approaches  infinity  for  either 
system,  this  limiting  value  is  approached. 

4.    OPTIMIZATION  STUDY  OF  COMPOSITE  SYSTEM 

From  the  standpoint  of  stability,  use  of  a  completely 
mixed  tank  reactor  followed  by  a  plug  flow  reactor  may  also 
be  desirable.   Because  of  this,  we  wish  to  investigate  the 
behavior  of  a  composite  system  under  optimum  conditions. 
The  schematic  diagram  of  the  composite  system  is  shown  in 
Figure  23- 


(1)   Formulation  of  Probl 


em 


The  problem  is  to  find  the  optimal  allocation  of  in- 
fluent and  the  optimal  distribution  of  volume  for  carrying 
cut  a  waste  treatment  process  in  a  composite  step  aeration 
system  composed  of  a  completely  stirred  tank  reactor  followed 
by  a  plug  flow  reactor.   The  allocation  of  influent  and  the 
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Fig.  23.  Schematic  diagram  of  composite 


87 


distribution  of  volume  are  to  be  chosen  so  that  the  total 
aeration  system  volume  is  minimized.   The  objective  function 
of  this  problem  is  therefore  the  same  as  in  previous  systems 
and  can  be  represented  as 

S  =  xf    .  (48) 

The  mathematical  equations  which  describe  this  system 
have  been  given  previously  as  both  systems  have  been  dis- 
cussed before.   The  performance  equations  for  the  mixing 
stages  are  Equations  (12),  (13),  and  (14).   The  perfor- 
mance equations  for  the  reaction  stages  are  Equations  (15), 
(17),  and  (18)  for  the  completely  stirred  tank  reactor  and 
Equations  (36),  (38),  and  (39)  for  the  plug  flow  reactor. 

This  optimization  problem  can  be  solved  by  the  discrete 
maximum  principle.   The  adjoint  equations  and  Hamiltonian 
function  for  n  =  1  are  the  same  as  those  for  the  completely 
stirred  tank  reactor  system  and  are  given  by  Equations  (19) 
through  (28).   Those  for  the  plug  flow  reactor  system  are 
given  by  Equations  (41)  through  (44). 

The  computational  scheme  for  obtaining  numerical  solu- 
tions which  was  discussed  previously  can  also  be  used  to 
solve  this  set  of  non-linear  algebraic  equations.   The  com- 
puter program  which  has  been  used  to  obtain  the  numerical 
solutions  is  given  in  Appendix  J> . 
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(2)   Results  and  Discussion 

In  this  numerical  study,  the  kinetic  constants  are 
assigned  the  following  values:   Y  =  0.5  mg/mg,  K  =  100 
;/l,  kfL   =0.1  hr   .   The  following  additional  conditions 
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are  also  specified:   R  =  O.L,    *2   =  150  mg/l ,  x^  = 

;/l,  x;  =  800  mg/l,  x|N  =  80  mg/l,  ^  0 ?n-l  =  4500  l/hr. 
*  n=l  1 

The  numerical  result  which  gives  the  optimal  influent  dis- 
tribution and  the  optimal  value  distribution  for  this 
system  is  given  as  follows: 

6  J  =  1810  l/hr, 

e2  =  1790  liters, 

ej  =  2690  l/hr, 

G^  =  6956.5  liters, 
V  =  e^  +  6^  =  8746.5  liters. 

This  result  may  be  compared  with  the  optimal  result 
fo  V  =  8358.1  for  two  plug  flow  reactors  connected  in 
series.   The  reason  that  the  result  obtained  using  the 
composite  system  is  not  as  good  as  that  obtained  using  the 
two  plug  flow  reactor  system  can  be  interpreted  in  Figure  24. 

•e  2L   the  optimal  locus  of  the  compostie  system  is 
represented  by  the  curve  ABCDEF.   If  for  the  same  operating 
conditions  the  first  completely  stirred  tank  is  replaced  by 
a  plug  flow  reactor,  the  optimal  locus  will  be  the  curve  ABDEF 
shown  in  Figure  23;  however,  this  is  not  the  optimal 
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.Fig.  24.  Growth  rate  trajectory  in  an    optimal    composite 
step    aeration   system. 


90 


condition  for  two  plug  flow  reactors  connected  in 
series.   As  described  previously  the  area  BCD  will  be  the 
hoi     Lirr.e  saved  by  using  two  plug  flow  reactors  rather 
than  using  the  composite  system. 

The  reason  for  often  employing  a  tank  reactor  as  the 
first  reactor  for  an  autocatalytic  reaction  has  been  dis- 
cussed by  Bischoff  (11).   It  is  used  in  an  effort  to  in- 
crease the  concentration  of  organisms  and  thus  to  attain 
the  maximum  reaction  rate  in  the  front  part  of  the  system. 
In  our  study  a  sufficient  amount  of  organisms  are  con- 
tained in  the  recycle  stream  and  thus  it  is  not  necessary 
to  use  a  tank  reactor. 
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CHAPTER  V 
CONCLUSION 

In  this  study  we  have  proposed  several  flow  systems 
for  use  with  the  step  aeration  process.   These  systems 
include  several  completely  stirred  tank  reactors  in  series, 
several  plug  flow  reactors  in  series,  and  a  composite 
system  composed  of  a  stirred  tank  reactor  and  a  plug  flow 
reactor.   Under  simplifying  assumptions,  the  optimum  oper- 
ating conditions  have  been  determined  by  using  the  maximum 
principle  and  a  search  technique. 

In  general,  from  the  optimization  studies  the  alloca- 
tion of  feed  along  the  system  is  found  to  be  desirable 
under  the  conditions  that  are  considered  in  this  study. 
In  other  words,  the  step  aeration  system  is  theoretically 
desirable  in  some  types  of  operating  situations. 

From  all  the  systems  we  have  considered,  we  have  found 
that  a  continuous  feed  allocation  along  a  plug  flow  system 
is  the  optimal  policy  for  the  step  aeration  system.   Under 
the  situations  and  assumptions  of  this  study,  the  plug  flow 
system  is  far  better  than  the  completely  mixed  flow  system. 
We  have  also  found  from  the  study  of  the  composite  system 
that  plug  flow  is  most  necessary  at  the  back  end  and  that 
complete  mixing  does  not  increase  the  volume  requirements 
very  much  at  the  front  end  of  the  step  aeration  system. 
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The  numerical  data  and  the  kinetic  constants  employed 
in  this  study  may  not  be  of  practical  significance.   However, 
. .::er  experimental  data  and  more  refined  kinetic  models 
can  be  easily  incorporated  into  the  scheme  developed  in 
this  study  to  obtain  more  realistic  results.   A  more  real- 
istic flow  model  can  also  be  used  instead  of  the  idealized 
reactors  we  have  considered. 
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NOMENCLATURE 


Symbol 

2 

A  Cross-section  area,   cm 

Hn  Hamiltonian  function  at  the  nth  stage. 

ks  Maximum  growth  rate  when  the  organic s  is  unlim- 

max  ited,   hr-1  . 

K  Growth  rate  constant,  mg/liter. 

L  Total  length  of  a  plug  flow  reactor,   cm  . 

£  Length  of  the  reactor,   cm  . 

r  Growth  rate,   mg/l-hr  . 

R  Recycle  ratio 

S  Objective  function 

s  concentration  of  organics,   mg/liter  . 

x  Concentration  of  organisms,   mg/liter  . 

x?  Volumetric  flow  rate  of  the  flow  discharged 
from  the  nth  stage,   l/hr  . 

xt>  Concentration  of  organics  in  the  inlet  stream 
of  the  nth  stage,  mg/liter  . 

x?  Accumulated  volume  up  to  the  nth  stage,   liter  . 

x^  Concentration  of  organisms  in  the  inlet  stream 

^  of  the  nth  stage,   mg/liter  . 

x*  Concentration  of  organics  in  the  feed,   mg/liter 

x.  Concentration  of  organisms  in  the  feed,  mg/liter 

Y  Yield  factor 

z?  The  ith  adjoint  variable  at  the  nth  stage. 

ex  Correction  factor 

&1  Volumetric  flow  rate  of  sewage  introduced  at  the 
nth  stage,   l/hr  . 

9 2  Reactor  volume  of  the  nth  stage,   liter  . 


PART  TWO 

AN  OPTIMIZATION  STUDY  OF  THE  REVERSE  OSMOSIS 
DESALINATION  PROCESS 
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CHAPTER  I 
INTRODUCTION 

Because  of  its  simplicity  and  the  potential  low  energy 
requirement,  reverse  osmosis  has  been  drawing  widespread 
attention.   Current  research  activities  have  been  directed 
to  (A)  development  of  membranes  which  combine  good  salt 
rejection,  high  water  transmission  at  reasonable  pressure, 
and  long  membrane  life,  (B)  improvement  of  membrane  fabri- 
cation techniques,  (C)  improvement  of  cell  design  and  con- 
struction, and  (D)  mathematical  and  experimental  studies 
of  salt  build-up  on  the  membrane  surface.   However,  only  a 
rather  limited  attention  has  been  directed  to  the  system 
analysis  of  the  entire  system  (12). 

Recently,  because  of  the  thermodynamic  undesirability 
of  mixing  a  recycle  stream  of  high  salt  concentration  with 
another  stream  of  much  lower  salt  concentration,  a  multi- 
stage sequential  system  with  brine  recycle  at  each  stage 
has  been  proposed  (13)-   In  order  to  increase  the  efficiency 
of  energy  recovery  of  this  system,  use  of  a  flow  work  ex- 
changer has  been  proposed  by  Cheng  et  al.  (14). 

The  present  study  is  directed  to  the  mathematical 
modeling  and  optimization  of  the  entire  reverse  osmosis 
system  incorporating  the  improvements  mentioned  above.   A 
combination  of  the  discrete  version  of  the  maximum  principle 
and  the  steepest  descent  search  technique  is  used  to  carry 
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it  the  optimization .   This  study  illustrates  how  optimi- 

n  car.  be  used  in  a  process  design  study  which  has  an 

objective  function  based  on  the  economics  of  the  process. 
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CHAPTER  II 
DESCRIPTION  OF  PROCESS 

1.    REVERSE  OSMOTIC  PROCESS  (15) 

uch  interest  has  been  shewn  recently  in  reverse  osmosis 
as  a  promising  method  for  the  economic  recovery  of  water  from 
saline  sources.   In  reverse  osmosis,  desalination  is  achieved 
by  forcing  salt  solution  under  pressure  through  a  membrane 
which  generally  passes  water  much  more  readily  than  salts. 
The  simplicity  of  the  reverse-osmosis  process  is  apparent 
upon  inspection  of  a  flow  sheet  of  the  Aerojet  pilot  plant 
(16).   Basically  the  process  consists  of  a  pumping  system  to 
raise  the  pressure  of  the  brine  solution  and  of  an  array 
of  selective  membranes.   The  only  energy  consumption  re- 
quired by  the  process  is  that  for  driving  the  pumps.   A 
reduction  in  energy  consumption  will  reduce  the  cost  of 
the  fresh  water  produced. 

Water  passing  through  the  membrane  is  supplied  to  the 
membrane  boundary  by  bulk  flow  of  solution  normal  to  the 
membrane.   Salt  is  carried  along  with  the  water.   If  a 
steady  state  is  to  be  maintained  without  an  accumulation 
of  salt  on  the  membrane,  this  salt  must  diffuse  back  into 
the  main  bulk  solution.   A  salt  concentration  gradient  is 
established  near  the  membrane  boundary  such  that  the  net 
salt  flux  normal  to  the  membrane  is  zero.   This  means  that 
the  salt  concentration  and  osmotic  pressure  are  greater 
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t  the     rane  thar.  in  the  bulk  solution.   If  we  use  a 
circulation  pump  to  obtain  fully  developed  turbulent  flow  ■ 

side  the  tubes,  the  bulk  of  the  brine  flowing  parallel 
tc  the  osrr.otic  membrane  surface  will  be  well  mixed,  and  the 
existence  of  concentration  and  velocity  gradients  will  be 
restricted  to  the  laminar  boundary  layer  as  shown  in  Figure 
1. 

Based  upon  this  boundary  layer  model,  the  desired 
performance  equation  which  relates  the  applied  pressure, 
circulation  rate,  and  the  flux  of  the  fresh  water  pro- 
duced has  been  developed  (13,  15)  in  detail. 

2.    DESCRIPTION  OF  MULTISTAGE  PROCESS  (13) 

As  mentioned  in  the  introduction,  because  of  the 
thermodynamic  undesirability  of  mixing  a  recycle  stream 
ci   high  salt  concentration  with  a  main  stream  of  lower 
salt  concentration,  a  multistage  process  has  been  pro- 
posed (13).   In  order  to  increase  the  efficiency  of  the 

rification  caused  by  the  stagewise  increase  in  salt 
concentration,  the  use  of  a  pressure  pump  between  suc- 
cessive stages  has  also  been  suggested.   Figure  2  illus- 
trates the  proposed  sequential  purification  process. 

The  process  is  a  form  of  a  simple  sequential  multi- 
st*    model  as  shown  in  Figure  1-13.   Each  stage,  except 
t  one,  consists  of  a  membrane  separator  unit,  a 


= 


99 


• 


^-  ^.  .  —  .    IU 


z' 


Ve 


- 


-      ,■-,- 


~/-/-/-/-/ 


boundary    Icysi 


. 


..     -   


. 


\ 


•^ 


V 


/ 


^v^..»y.  >-/•  »_•_.. 


-  -  •-  , 


-    -       .     -^  -        •--         n        ---  —         _   -    -    . 


: 


c 

v. 
- 

O 

a 

(/) 

r « 

> 
o 

CJ 

' 

Z  cvj 

1 

-  ?3 

C   OJ 

4 

C/3 

c 
L 

■_ 

'.' 


<v  CM 
Til     r^" 


—  rv 


CV 


til  A 


CJ 

3 


-  u 


UJ  cf 


M 


r> 


— 

,2  - 


»~^ 

(    •        r— 


rj 


C 


_ 


ri  /-• 


' 


- 


OJ 


101 


recycle  unit  and  a  high  pressure  pump  unit.   The  last  stage 
is  similar  to  the  others  but  adds  one  turbine  unit  in  its 
outlet  for  energy  recovery.   When  a  flow  work  exchanger 
is  incorporated  into  the  proposed  sequential  process;  it 
is  connected  between  the  first  and  the  last  stage  as  shown 
in  Figure  3 • 

In  setting  up  these  models,  the  following  assumptions 
have  been  made. 

(a)  The  flow  model  is  fixed,  but  the  geometric  par- 
ameters of  the  membrane  separator  unit  itself,  i.e.,  the 
pipe  diameter  d,  length-to-diameter  ratio  L/D ,  and  the 
total  number  of  tubes  used,  m,  are  chosen  arbitrarily. 

(b)  The  salt  concentration  in  the  fresh  water  pro- 
duced is  assumed  to  be  zero,  i.e.,  the  salt  cannot  pass 
through  the  membrane . 

(c)  The  saline  water  feed  solution  contains  3«5 
weight  per  cent  salt. 

(d)  There  is  no  precipitation  of  salt  on  the  membrane 
surface. 

(e)  The  costs  of  the  pump,  turbine,  and  motor  are 
assumed  to  be  directly  proportional  to  horsepower  rating 
in  the  horsepower  range  of  interest. 

(f )  The  cost  of  the  membrane  separator  unit  is  pro- 
portional to  the  weight  of  the  material  used. 

(g)  Da ,  jj.  ,  andpare  assumed  to  be  constant  in  the 
concentration  range  of  interest. 
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(h)   Membrane  constant  K  is  assumed  to  be  independent 
of  pressure. 

Process  and  cost  analyses  for  the  multistage  reverse 
osmosis  process  have  been  developed  previously  (13).   We 
now  simply  summarize  these  relations  as  follows: 
(I)   Process  Analysis 

The  equation  for  the  fresh  water  flux  at  the  nth  stage 
is 

Pn  .  KUpn  -  12100  xn] 

—   —  

_n 
1  +  c   x 

(Ren)7/8 
where 


c  =  3.05x10^  -573 


(Sc)"Da 


K  =  the  membrane  constant  — ^  - 

ft  -hr-psi 

a p1^  the  pressure  difference  across  the  membrane  of 

the  nth  stage,  psi. 

Ren  =  the  average  Reynolds  number  at  the  nth  stage, 

S  c  =  S  c  hmi d t  numb  e  r , 

d  =  the  diameter  of  the  tubes  in  the  membrane  separator 

unit ,  ft , 

Da  ■  molecular  diffusion  coefficient  of  salt,  Sg  c~  , 
a  '   bee 


x"  ■=  the  average  mass  fraction  of  salt  at  the  nth  stage. 

z:'f.e   overall  material  balance  and  the  salt  balance 
around  the  nth  stage  give 

,»n  -  „n-l   v-n  .  „n-l   vri  cn  /0\ 

q-q     -..   -q     -  1*  pS   ,  [di 

c°x°  =-  qnxn  ,  (3) 

where 

c/'>  qn~^  =  the  mass  flow  rate  of  the  brine  solution 

discharged  from  the  nth,  (n-1)  stage,  lbn. 

hr 

q°  =  the  mass  flow  rate  of  the  inlet  brine  solu- 

lb„ 

tion,  S  . 

hr 

Wn  =  the  mass  flow  rate  of  fresh  water  product 

1  "hi 

from  the  nth  stage ,   m  . 


hr 

lb. 


P  -  the  density  of  the  fresh  water, 


ft3 

sn  ■  the  membrane  area  of  the  nth  stage,  ft2. 
x°  -  the  inlet  mass  fraction  of  sea  water. 

Combining  Equations  (8)  and  (9)  and  using  a  salt  balance  to 
eliminate  qn~l  gives 

.ovo       vo  „n-l 


x     x  (/J 

x"-xi4  *  F"   ,Sr 


.n     °   n-±  „n   IH 


p(7o) 


q< 
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(II)  Cost  Analysis 

The  total  cost  required  for  the  system  mainly  consists 
of  the  operating  and  capital  costs.  These  costs  are  separ- 
ately summarized  in  Section  A  and  B  below. 

A.   Operating  cost: 

The  energy  requirement  for  the  high  pressure  pump  at  the 
nth  stage  per  pound  of  total  fresh  water  produced  is 


E-/  1  t  0  «   n   .  n-1         o 
'  f  ap  -  S  p      x~ 


"f  <7m  Up      P       xn-l  (1.  4, 


(5) 


xK 


where 


*?  ,7  ,7  =  the  mechanical,  pump,  and  turbine  efficien- 
m  p  r  7.1-j 


cies , 


™f   =  the  loss  factor, 

Wf  =  Total  mass  flow  rate  of  fresh  water  produced 
from  the  system,  lbm/hr  . 

The  energy  requirement  for  the  recycle  pump  at  the  nth 
stage  per  pound  of  total  fresh  water  produced  is 

E2   =  0  023  (Re*)2'8  (i£)3      SP       1  +  ^f 
d        U.O^  [rte    )  (       )  (o) 

Wf  p    qo  (i  _  *)    ^c'm'p 

The  energy  recovery  at  the  reject-brine  turbine  per 
pound  of  total  fresh  water  produced  is 
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-  =  V7n  U  -  ^J  ^    VX°   o  (7) 

p   m  f     e     Xs  -  x°  • 


If  a  flow-work  exchanger  is  used,  the  turbine  is  only- 


used  to  recover  energy  as  the  pressure  is  reduced  frorr.  Ap 


:■: 


to  Ap  .   The  flow  work  exchanger  is  in  turn  used  to  transfer 
energy  to  the  feed.   An  amount  q^   of  sea  water  is  pressurized 
ir.  the  first  stage  with  this  energy.   Therefore,  the  energy 
recovery  per  pound  of  total  fresh  water  produced  becor.es 


E->  \^T  -  A^1   v° 

r^  =  V     <?   (1  -  I?  )  — ±2 2£ (g) 


Because  of  the  flow  work  exchanger,  the  energy  requirement 
for  the  high  pressure  pump  at  the  first  stage  is  simply  to 
pump  the  amount  of  q°  -  c  rather  than  q°.   Therefore, 


1  _               1      An 
"f        ?rr.  ^p          e 

0          N 

q    -  q 

o.r.ce 

Wf   =    qO    _    qN    , 

Apc    =    0, 

(9) 


Equation  (9)  becomes 


^r^%  ~  ■  (10) 
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Here  all  the  energy  required  is  assumed  to  be  supplied 
from  electricity  and  the  electrical  power  cose,  C  ,  is  as- 
sumed to  be  $0,007  per  Kw-hr  for  all  cases. 
B.    Capital  cost: 

For  simplicity  the  cost  of  pump,  turbine,  and 
motor  are  assumed  to  be  directly  proportional  to  horse- 
power rating  in  the  horsepower  range  of  interest. 

The  membrane  separator  cost  is  assumed  to  be  pro- 
portional to  the  weight  of  the  material  used.   The  weight  of 
the  material  used  per  pound  of  the  fresh  water  produced  for 
this  unit  has  been  suggested  to  be 

4t^  ^  +  ^+2#2/5)    (ii) 


wnere 


?      -   the  density  of  the  material  of  construction,  ^m  , 
m  3 

ft 

Wg  =  the  mass  of  the  shell-and-tube  membrane  separator 
unit  of  the  nth  stage,  lbm, 

°~is  =  the  allowable  stress  of  the  material  of  construction, 
psi, 

L/D  =  the  overall  length-to-diameter  ratio  of  the  mem- 
brane separator  unit. 
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In  the  final  cost  estimation  by  Lonsdale  et  al.  (15), 
.b.  costs  of  $50/Kw  were  assumed  for  the  circulation  pump 
h  motor,  ar.d  the  blowdown  turbine,  an  f.o.b.  cost  of 

..  was  assumed  for  the  high-pressure  pump,  and  an  f.o.b 
cost  of  s'4.40/lb  for  the  separator  unit.   As  recommended  in 

-  Dffice  of  Saline  Water  Procedures  (17),  f.o.b.  costs 
should  be  multiplied  by  1.3  to  obtain  installed  costs.   It 
also  recommended  that  interest  on  the  depreciated  plant 
investment  be  at  4%,  and  depreciation  be  over  20  years. 
This  gives  a  total  charge  of  7.4%  for  amortization  on  the 
initial  investment.   In  addition  to  the  above  items,  two 
per  cent  is  added  for  taxes  and  insurance.   This  gives  an 
annual  capitalization  charge  for  these  equipment  items  of 
9-U%   of  the  initial  cost.   Then  an  assumption  of  a  load 

factor  of  330-on-stream  days  per  year  gives  a  capitalization 

-6 
charge  of  11.94x10   of  initial  cost  per  hour  on  stream,  ^. 

C.   Total  cost: 

Combining  both  the  operating  and  capital  costs, 

the  total  cost  contribution  of  the  system  is  obtained  in  the 

form 

rr  En  N   E^ 

CT  «  <  <V  +  Ce)  £  £   +   (VCpR  +  Ce)  £   ^ 


+  (yc  -  c  )  I2L  +<pc  y       s  <12) 

pT    e  Wf   YUs  ^  W7 
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whore  C  ,.,  C   .Cm,  and  C   are  respectively  the  capital 
pH'   pR   pT"      s 

costs  of  the  high-pressure  pump,  circulation  pump,  turbine, 

and  separator  unit. 

Substituting  Equations  (5),  (6),  (7),  and  (11)  into 

Equation  (12)  and  combining  all  constants,  we  obtain 

o    N   .  n  A    n-1       N 

1  -  Tf  xx 

x 


r  r  2  3   sn        N    n  sn        N     n  *  sn  1 

•B2  SI  (Re11)2*8  (%)  +  B3  Z  APr'(%)  +  B4  r  Upn)  (%) 

^   n=l  4        n=l     q        n=l       4 


x°    _N 
H  -  x( 


+  B5^V^4p  (13) 


where 

1  ?H    e    f?m^p 

B9  =  0.023  (YC  -  +  Cj   — - — £   (4r)3 

2  pR    6   ScVp   <*? 

B3 (1.62  +  -j-^)    , 

m  f 

B  -°'189  W 

4    ffm  L/° 


35  =  «iC  T  -  Ce)  i?p  ^m  (1  -  f?f)/?  . 
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hen  a  flow  work  exchanger  is  employed,  we  can  combine 
ions  (6),  (S),  (10),  (11),  and  (12)  to  obtain 


B-,  -  ?" 


o    N 

1  ,    2^  n=2 
1  "  N 


Apn  -  A  P 


n-1 


n-1 


(--JE )   B0  r  (Ren)2'8(2_)-fB,  r    ^pn(^-) 

Xs  -  x°  I  2  n-1         q°   3  n-1     q° 


r 

n=l 


i.-n.        xo 


B.  H   Upn)$(^r)   +  B^ 

q°  J     >  x\-  x° 


UpN  -  Ap1) 


-  3,  -*_  pN   . 

o  ^r  xo 


(14) 


where 


O 


_  1 


£  WCpF+SCeJ 


re  we  assume  that  the  capital  cost  of  the  flow  work  ex- 
changer, C  p,  is  the  same  as  that  of  turbine  and  use 

to  denote  the  inefficiency  of  the  device. 
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CHAPTER  III 
OPTIMIZATION  STUDY  OF  THE  PROCESS 

In  Chapter  II  the  equations  which  describe  the  perfor- 
mance of  the  reverse  osmosis  process  are  given.   Cost  equa- 
tions which  relate  these  performance  equations  to  the  unit 
cost  of  product  water  are  also  given. 

In  this  chapter  these  performance  equations  and  cost 
equations  are  used  together  with  optimization  techniques 
to  determine  the  values  of  the  design  variables  that  mini- 
mize the  cost  of  producing  a  given  quantity  of  fresh  water. 

In  the  present  study  we  classify  the  optimization 
studies  of  the  reverse  osmosis  system  into  the  following  two 
classes  of  problems: 

(1)  Multistage  operation  without  the  use  of  a  flow  work  ex- 
changer, 

(2)  Multistage  operation  with  the  use  of  a  flow  work  ex- 
changer. 

As  in  Part  I,  the  discrete  maximum  principle  is  also 
employed  here.   The  discrete  maximum  principle  has  been  briefly 
described  in  Part  I. 

1.    MULTISTAGE  OPERATION  WITHOUT  FLOW  WORK  EXCHANGER 
(1)   Formulation  of  the  Problem 

The  optimization  problem  will  be  formulated  so  that  it 
conforms  to  the  "standard"  notation  and  structure  of  a  dis- 
crete maximum  principle  optimization  problem.   Figure  4 
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illustrates  the  structure  of  the  sequential  three  stage 
discrete  maximum  principle  model  of  the  reverse  osmosis 
system  without  the  flow  work  exchange  that  is  shown  in 
Figure  2.   In  Figure  4  two  types  of  variables  are  indi- 
cated.  Those  represented  by  are  decision  variables  which 
are  treated  as  independent  variables,  while  those  repre- 
sented by  x  are  state  variables  which  are  dependent  vari- 
ables.  Superscripts  are  used  to  denote  the  stage  with 
which  each  variable  is  associated.   The  state  and  decision 
variables  are  defined  as  follows: 

(a)  State  variables 

x?   =  salt  concentration  at  the  nth  stage, 

x  =  total  cost  required  up  to  and  including  the 
nth  stage. 

(b)  Decision  variables 

n     n 
G]_  =  Apl  ,  the  pressure  difference  across  the  mem- 
brane of  the 'nth  stage, 
©2  =  Hen,  the  average  Reynolds  number  at  the  nth 

stage, 

sn 
0  ?  =  —p;   >    the  ratio  of  membrane  area  and  inlet  mass 

->   q 

flow  rate  at  the  nth  stage. 

From  Equation  (2.4)  we  can  see  that  the  performance 
equation  for  the  first  state  variable,  x-j_,  can  be  expressed 
as 


-i; 


...  _  ....  j. 


3o 


yn 


-  12100 


n 


n  -  1,  2, 


(eS) 


7/8 


(D 


„re 


yn  = 


.n 

'1 


(2) 


The  performance  equation  for   the  second  state  variable 
x'p",  can  be  obtained  from  Equation  (2.13)  as 


r.    n-1 


n-1  //.n   „  n-1 


+  A-j_y    (6^  - 


n 


)  +  A26o 


B2(G2)2*8  -  B3(eJ) 


+  E 


4 


(ej)*)  , 


n  =  1,  2,  . . . ,  N-1, 


(3) 


:;   ::-i   .  n-i  ,.n   .  ::-i.   a  if  ,flN.2.8   n  /flN. 

x2  =  x2   +  AQ_y    (t?i  -  0]_   )  +  A2^3   B2(82)    +  33(01) 


. 


+  b4(0JJ)2     +  A3ej   , 


N 


(4) 


Since  Equations  (3)  and  (4)  have  memory  in  decision 


,  for  simplification,  we  can  introduce  the  new  state 
r. 

variable,  x?  such  that 


n   ,  r 

-  e     , 


'J 


n  =  1 ,  2 ,  .  .  .  ,  N 


(5) 


aations  (3)  and  (4)  then  become 
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n-1   »   n-1  /nn    n-l\    .  Qn  ~  /,n\2.8   n^n 
x0  ■*  x«   +  A-j^y    [S^   -  x2      '   +  A2C3   32^2^       3^1 


. 


+  3^(6^   ,      n  =  1,  2,  ...,  N-1,  (6) 


::    N-1    .   N-1  ,„N    N-1,    .  AN   „  ,nN,2.3    „  aN 


x2   x2 


iy^-x  (ej  -  x?    )  +  A2e3    b2(©2)       +  B3e1 


+  34(6^)i    +  A3eJ 


where 


(7) 


Al  ~  Bl  ZG  (3) 

xl 


-1 

A2  "  ~N "0  (9) 

Xl  "  Xl 


A3  =  33?^?   •  (10) 

xl  "  xl 


are  to  find  the  optimal  design  variables  that  mini- 
mize the  total  required  cost  of  producing  a  given  quantity 
of  fresh  water.   The  objective  function  of  this  problem  is 
therefore  simply  represented  as 

S  =  x|T  .  (11) 
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ying  a  discrete  from  of  the  maximum  principle 
:c  determine  the  optimal  policy  for  the  system  as  described 
above,  the  following  Hamiltonian  functions  and  adjoint  are 
roduced 


=  z£yn  +  z£  xlj  +  zn  xn  ,  (12) 

n-1    __1    n  -z>  yn     n*x2      n  3  A      ,  (1^) 
Zl    ay^"1   Zl  ay^-1  "  Z2  »yP-l  T  z3  ayn-i 

x2       aX2       aX2      ox2 

n         n 

n-1  _  3  Hn  .  zn  a  yn  +  zn  5X  2  +  zn  ^  x3   .  (15) 
3    ax-"1    1  Sxn-1    2  ^n-1    3  Bxn-1 

Differentiating  Equations  (1),  (6),  (7),  and  (5)  with 

respect  to  y11"1,  x§~  ,  and  x!?~  ,  and  substituting  the  results 
into  Equations  (13),  (14)  and  (15)  respectively,  we  obtain 

zn-l  .  zn  (G)2 +  7n  ,  /fin  Yn-!) 

21     "1              ~ +  z2  V91  X3   j 


2  +  kie^T7nT77^+12,100] 


(or  +  a  ^  (6?) 

2  (16) 


jre 


G  "  ^  +   „%;-  '  n  -  1,  2,  ...,  N,  (17) 
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n-1 


+  z?  ,   n  -  1,  2,  . . .,  N, 


(18) 


n-1  =  ,n.   n-1 
-3     -2Aly 


n  =  1 ,  2 ,  •  .  •  ,  N 


(19) 


N 
Since  x^  is  specified,  the  boundary  conditions  for  the 

adjoint  variables  in  this  system  are 


zN  -  1 
Z2    ±} 


(20) 


N   n 

Z  .}   —  U. 


(21) 


According  to  the  so-called  weak  form  of  the  discrete  maximum 
principle,  the  optimal  decision  variables  at  each  stage  must 
satisfy  the  equations 


->  rTn 

-.    n                 3x5             sx^ 

zn  dv        +  2n          2   +   2n  _1  = 

1  ae  n       2  dbn       3aen 
i               i              i 

o, 

i 
n 

■ 

1,  2,  3-     (22) 

-i-  j    *C  y         m   •  •  |    i  »  • 


where 


-K,y 


n 


ae£ 


vn  +      c    i  nn 


y? 


(eg) 


£77^  I   +  Kl  0  3 


cG 


(e?)7/8 


+  12100 


n-1,  2,  .  .  .  ,  N. 


(23) 
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I   v 


n 


-  n    /..n  ^n 


-7  ck  1  e?  (yx  ex  -  12100) 


n 


•     L5/8  f /  n  c 

J  (yx  -  _ 

n,7/o 


n    n 


K1e3e1  (yj 


(etp 


(«2) 


.) 


.n-  n 


-K^^yVC-V  -   12100)1 


n  =  1 ,    2  ,    .  .  .  ,    N , 


(24) 


,n      n 


Kx    O"  y     -   12100) 


r. 


yn  +       c 


+Kle3 


n 

c  1 


+  12100 


.n 


TifT^ 


(25) 


1*  JL   «         *  -    >  •     •     •     j  -  *    • 


a^11"1  +  a2b3  e^  i  a2b4  e^(ej)*  , 


n  =  1,    2, 


K  -   1 


(26) 


se 


Al,/"""L  +  A2B3  e3+  i  A2B^(qI)"*  +  A3    , 


(27) 


r. 


2.8  A2B2    (0")    (9g)' 


1.8 


n  ■  1 ,   2 ,    . . . ,   N , 


(28) 


lie 


-^  -  A2  i32(e5)2*B  T  B^  e?+  3,  W?)*l  .   n  =  1,  2,  ...,  N. 
-^ « n    *  J.  *   2        3   1   4  X   J 

3  (29) 


axn 


36? 


2.-1,  n  -  1,  2,  ...,  N.  (30) 


?  n 

_2  =  0,  n  -  1,  2,  ...,  N.  (3D 

,  s  n 

c-o2 


Bx11 


1_  =  0,  n  =  1,  2,  ...,  N.  (32) 
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(2)   Results  and  Discussion 

The  computational  scheme  for  carrying  out  the  optimi- 
zation of  the  sequential  multistage  problem  with  some  end 
points  fixed  as  described  in  Part  I  can  be  applied  to  this 
problem. 

Numerical  computation  has  been  performed  with  an  IBM 
162C  computer.   A  computer  diagram  to  carry  out  a  three 
stage  problem  is  shown  in  Figure  5.   A  computer  program  is 
given  in  Appendix  4- 

The  following  data  are  chosen  for  this  study: 


(start  \ 


v3     .      ©'       01       ft' 


ef  /  e| ,  e|  ,  e|  , 


:     ' 


ER  S;  AQn    and   cc 


© 


VI 


y2     vn 


,  c 


t      >      ftO 


it 


e?  -  en-  A9n 

2  •      3  '      I   »     t» 


aej, 


Qj-e^+AeJ 


-':-.. j.    Cc-^1.";-:'      *:?•  ; 
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-5   2 

D     =  1.5  x  10   cr"-   ,  the  diffusion  constant  for 

sec 

sodium  chloride  in  water  at  room  temperature, 

K     =  0.86  x  10_z*-  — pi , 

ft2-hr-psi 

1)     =  o  x  10" ^  cm£  for  the  brines  at  30°C, 

sec 


p     =  62 -4 


'El 


lb"B 


ft3 

=  i   f  t .  , 
24 


D     ~  8> 


lb 

=  2.7  x  62.4  — 7  ,  density  of  aluminum, 
ft-* 


0"  -  15000  psi,  allov/able  stress  of  aluminum, 

^  =  11.94  x  10"°  hr"1, 

C  „  =  100  -f-  f .o.b., 

pH  KVJ 

C^T  =  50  -^  f .o.b. , 

?T  KW 

CnR  =  50  -&-  f .o.b., 

pR  '  KW 

Cg  =4.4  yJc —  of  aluminum  f.o.b., 

^Dm 

Cp  =  0.007  — $ —  , 

e  Kw-hr 

?n  =  0.9  , 

%  =  ^?p  =  0.8  , 

f?f  =  0.1. 
=  0.1. 
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Tie  -     ^  constants  are,  therefore, 


1 

1 

-  (iL)3  -  8.1, 
Da 

0 

«  >-0^  1Q5  Kd  =  2320, 
Sc1/3D  a 

Bl 

=  0.1136  x  10"7   $ 

psi-lbm 

32 

=  0.6330  x  10"17   $ 

ft2-hr 

33 

=  0.52  x  10"7  — | —   , 
ft-  -hr 

3 

=  0  177^  y  in- 6      $ 

4 

psi-ft-hr 

B5 

=  -0.3494  y  10"8   $ 

lb  -psi 

36 

=  0.11238  y  10"8    <c: 

lbm-?si 


The  numerical  results  are  listed  in  Table  1.   The  op- 
1  cost,  optimal  Reynolds  number,  optimal  pressure  drop, 
imal  membrane  area  to  feed  rate  ratio  are  plotted 
;ainst  the  parameter,  y^,  in  Figures  6,  7,  8,  and  9-   From 
■chese  figures  we  can  find  the  overall  optimal  cost  and  the 
ccrresp      5  design  variables.   They  are 

•  -  U01.4,    A?2  -  1686.9,     ^P3  =  1704.3, 
,2.0,    Re2  =  9088.0,     Re3  -  10030.0, 
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1 


2  c3 


~-  0.0533,     ^=0.0524,     S_=  0.0438, 

q  qu  q 

=  0.09 

Total  cost  =  CO. 43 97/1000  gal. 

r  one  thousand  gallons  of  fresh  water  produced  per  hour, 
the  feed  rate,  q  ,  can  be  calculated  from  the  overall  ma- 
terial balance  and  the  corresponding  membrane  areas  of 
individual  stages  can  be  evaluated.   They  are 

q°  =  13619.5  lbm/hr., 

s1  =  725.92  ft2,   s2  =  713.66  ft2,   s3  =  596.53  ft2 

2.     ULTISTAGE  OPERATION  WITH  FLOW  WORK  EXCHANGER 

(1)   Formulation  of  Problem 

The  problem  considered  here  is  an  extension  of  the 
pro!     in  the  preceeding  section.   All  the  performance 

-ions  are  the  same  as  those  for  the  previous  problem 
except  for  the  cost  functions  as  analyzed  in  the  previous 
.  .   Therefore,  in  this  section  we  only  need  to  refor- 
Late  the  performance  equation  for  the  second  state  vari- 
The  objective  function  is  again 


(33) 
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From  Equation  (2,55)  the  performance  equations  for  the 
second  state  variable,  x£,  can  be  expressed  as  follows: 

x2  =  x2  +  A2  e3  [B2(e2)2*8+  B3  Gl  +  h{el]    ]»      (34) 


n 


n-1        .        n-1    /nn         n-l»    ,,     fln  fD    ,Qn» 
=  x9        +  A-l  y  ^61  "   x3      '      A2     3   I  22  ^    2 ' 

+  BLeJ  +  B4(0j)*J    ,      n  -  1,   2,    ...,   N  -1, 


n ,  2 . 3 


(35) 


N  N-1        ,      N-1    /QN       fl N-1,         .    Jfn    /flN*2.8 

x2  =  x2  ly  ^    1  1      '  2  3      V    2' 


N 


N 


T.n 


+   33(G£)    +   B^    (e£)a}+        (A3+  A4)0i^    , 


>N 


(36) 


where 


A    ,    A    ,    and  A~   are   the   same  as   defined  before, 


h  =  B6 


X 


0 


u 


x- 


(37) 


Similarly,  the  Hamiltonian  functions  and  adjoint  vari- 
ables as  described  in  Equations  (12),  (13),  (1A-),  and  (15), 
and  the  necessary  condition  expressed  by  Equation  (22)  are 
valid  for  the  problem. 

The  derivative  terms  relative  to  the  second  state  vari- 
Le,  x^",  can  be  summarized  as  below: 
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dx 


3y 


n 


r.-l 


n-1 


=  A1(6J  -  x^"1) 


n  =   2,    .. . ,    N.         (33) 


Sx 


n 

2       =1 


n  =   1 ,    2 ,    .  .  .  ,    N , 


(39) 


dx 


n 


3  xV 


A  n-1 


n  -   1 ,    2 ,    .  .  .  ,    N . 


(40) 


3x; 


9y 


o 


-  o, 


(41) 


2x 


26 


-  A 


-LY11"1  -  A2B3G3   +   -iA2B46^    (e£)      ,      n  =   1,    2,    (42) 

•  •  •  a      IN  *  -»~ 


^  _ 


l      1 


ae? 


-  B1  +  A2B303   +   2  A2   B4S3    ^1*        "A 


(43) 


a  x 


ae 


A^N-1  +  A2   B36^  +  |  A2   B4ef(ejH   -    (A3   +  AJ 

(44) 


•  '   E2     V3(e2}1'8 


n-1,   2,    ...,   N.  (45) 


BQ 


n 


=  A2  [b2(9^)2-8  +  B3  ©  ^  +  B4(e£)r|   » 


n  =  1 ,  2  ,  .  .  .  ,  N . 


(46) 


(2)   Results  and  Discussion 

By  using  the  sane  computation  scheme  as  employed  in  the 
preceeding  section,  a  numerical  computation  was  performed 
with  an  IBM  computer  1620.   The  computer  program  for  his 
study  is  given  in  Appendix  5.   The  same  numerical  date  were 
also  used. 

The  numerical  values  of  optimal  results  are  listed  in 
le  2.   In  Figures  10,  11,  12,  and  13,  the  optimal  cost, 
optimal  pressure  drop,  optimal  Reynolds  number,  and  optimal 

membrane  area  to  feed  rate  ratio  are  plotted  against  the 

N 
parameter,  x-^.   From  these  Figures  we  can  fing  the  overall 

optimal  cost  and  the  corresponding  design  variables.   They 

are 


z.?1  =  1355, 
Re1.  =  6070, 


A?     =   1390, 
Re2  =  SS20, 


s   =  0.0465,    s£  =  0.0465, 


0 


A?3  =  1375, 
Re^  =  9470, 

s3  =  0.0465, 

q 


=  0.0065, 


Total  cost  =  $0.3636/1000  gal. 
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For  one  thousand  gallon  of  fresh  water  produced  p  . 
hour,  the  feed  rate,  q  ,  can  be  calculated  from  overall 
material  balance  and  the  corresponding  membrane  areas  of 
individual  stages  can  be  evaluated.   They  are 

q°  =  18,076  lbm/hr;    S1  -  S2  =  S^  -  840.53  ft2. 


By  comparing  the  optimal  results  for  the  processes 
with  and  without  the  use  of  the  flow  work  exchanger,  the 
following  observations  can  be  made. 

A.  The  use  of  the  flow  work  exchanger  gives  rise  to 
savings  of  16  per  cent  in  the  total  cost. 

3.  The  use  of  the  flow  work  exchanger  substantially 
reduces  the  energy  requirement. 

C.  The  outlet  salt  concentration  shifts  from  0.09 
mg/l  to  0.065  ^gA  • 

D.  Because  of  the  reduction  of  the  outlet  concentration 
the  feed  rate  increases  for  producing  the  fixed  amount  of 
fresh  water  when  the  flow  work  exchanger  is  used. 

E.  The  use  of  the  flow  work  exchanger  reduces  the  re- 
cycle ratio  at  each  stage. 

F.  The  Reynolds  number  for  flow  through  each  stage  de- 
creases only  slightly  when  the  flow  work  exchanger  is  used. 
This  is  due  to  the  two  compensating  effects  given  in  D  and  E. 

G.  The  use  of  the  flow  work  exchanger  increases  the 
membrane  areas  required  at  individual  stages.   It  is  also 
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.XT 

designs  for  a  step  aeration  activated  sludge  process  and 
osr.-.osis  desalination  process  are  carried  out  using  a  discrete 
I  principle  together  with  a  steepest  descent  search 
."our  types  of  systems  are  investigated  for  the  overall  pattern 
of    .  -ng,  and  distribution  of  feed  for  the  waste  treatment  process. 

are  developed,  which  enable  the  feed  allocation  and  the  holding 
to  be  chosen  such  that  the  total  required  volume  is 
use  of  a  plug  flow  reactor  and  the  allocation  of  feed  along 
..-e  found  to  be  theoretically  desirable  for  the  process.   Con- 
.ous  feed  allocation  along  the  plug  flow  reactor  is  found  to  be  the  best 
theoretically  for  carrying  out  the  activated  sludge  waste  treatment 
process.   Relations  for  the  reverse  osmosis  desalination  process  are 

..able  the  pressure  difference  across  the  membrane,  recycle 
ratio,  and  the  membrane  area  for  each  stage  to  be  determined  such  that  the 
sh  water  production  cost  is  minimized.   A  flow  work  exchanger  is' 
loyed  for  the  desalination  process  in  order  to  increase  the  energy  re- 
covery at)      .  reduce  the  total  production  cost.   From  the  analysis,  the 
.  ie  use  of  the  flow  work  exchanger  in  the  process  is  found  to  be 
. cant. 


